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ABSTRACT 

This thesis presents an analysis of transonic potential 
flow through an oscillating unstaggered thin plate cascade. 
A collocation technique is used involving the superposition 
of adjacent blade "isolated foil" potentials with inter- 
ference potentials of unknown strength. Imposition of 
flow tangency requirements leads to integral equations for 
the unknown source distributions of the interference 
potentials. Results presented include the interference 
coefficients for the steady and oscillating cases. 

The steady case is extended to the determination of the 
potentials along the blade and compares favorably with a 
parallel solution using a Laplace transform approach to 
the sonic wind tunnel wall interference problem. 
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I . INTRODUCTION 



Considerable advances are still being made in the 
development of aircraft turbopropulsion systems. In the 
quest for higher thrust to weight ratios, one must look to 
higher stage loading as well as significantly increased 
turbine inlet temperatures. From the compressor and fan 
design point of view, this translates primarily into increased 
tip speeds. Moving toward better blade design, the inves- 
tigation of transonic and supersonic relative flows is 
becoming more important not only for increased performance 
but also for flutter prediction. As shown by Pratt and 
Whitney, various flutter problems have already been encoun- 
tered [Fig. 1 as extracted from Ref. 1] and more such 
problems can be expected in the future. 

Recognizing that the actual three dimensional flow 
fields encountered in a transonic stage are complicated by 
the effects of rotation, finite blade thicknesses and cambers, 
boundary layers, and spanwise Mach distributions, it is 
excessively difficult to model these effects simultaneously. 
One can hope that eventually a complete analysis will 
become practical, but at this time the problems must be 
faced piecemeal. One simplified technique historically 
used in modeling such flows is the two dimensional cascade. 

The following investigation was concerned with transonic 
cascade flow. This problem had previously been considered 
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Figure 1. SCHEMATIC COMPRESSOR MAP SHOWING FLUTTER BOUNDARIES FOR FOUR TYPES OF SIMPLE MODE FLUTTER 



by Haraamoto [Ref. 2] and Gorelov [Ref. 3]. Also of impor- 
tance is the related problem of transonic wind tunnel 
interference. Here the studies of Sandeman [Ref. 4] and 
Murasaki [Ref. 5] must be mentioned. For a more general 
exposition of this problem, reference is made to Berndt 
[Ref. 6] and Goethert [Ref. 7]. 

Hamamoto based his analysis on the linearized oscilla- 
tory transonic flow equations. Gorelov attempted to include 
the effect of blade thickness by approximating the nonlinear 
transonic perturbation equation. This approximation, 
however, is valid only for low supersonic f lows . In the present 
investigation, a transonic flow approximation was intro- 
duced based on previous work by Oswatitch [Ref. 8] and Teipel 
[Ref. 9]. The technique permitted a general formulation 
for both the steady choke flow problem (in conjunction with 
Spreiter’s local linearization technique [Ref. 10]) and 
the oscillatory problem. The resulting boundary value 
problem was solved in a manner similar to that used by 
Gorelov. The problem formulation, solution technique, 
numerical evaluation of the recurring integrals, and 
representative results are described in the following 
sections. 
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II. UNSTEADY TRANSONIC FLOW THEORY 



A. DEVELOPMENT OF THE GOVERNING TRANSONIC FLOW EQUATIONS 






Y 



Uo» 







Figure 2. COORDINATE SYSTEM 



Consider planar flow of an inviscid adiabatic perfect 
gas along the positive X axis which is perturbed by a thin 
body (or thin bodies) executing small transverse unsteady 
motions. The perfect gas equation is applicable 

p’= pRT (II-l) 

where : 

p = pressure 
p = density 

T = absolute temperature 
R = gas constant 
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Furthermore the conservation conditions for mass, momentum, 
and energy must hold. Denoting the velocity components in 
the X and Y direction respectively as U and V, the conserva- 
tion of mass may be expressed as 



i£> + 3 (Up) + 9 (VP) = o 

at ax a y 



(H-2) 



Conservation of momentum is expressed 



du a_u , TT a_u , v au = _ i ap 

Dt at ax aY p ax 



(II-3) 



dv _ av TI av V av = _ l 3P 

Dt at ax aY p a y 



Conservation of energy is 



Dr . W n 
Dt + 2 ^ 



- p Cfx (up) + Iy ( vp)] 



(11-4) 



or 



Ds = o 
Dt u 



(II-5) 



where 



D_ 

Dt 



= substantial derivative 



e 

s 



= energy for the individual particle 



= entropy for the individual particle 



W = 



2 2 
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These equations, of course, exclude finite discontinuities 
(shocks) in the flow field. 

As is well known, an important role is played by certain 
functions which reduce the number of dependent variables. 

Most often the velocity potential is introduced. Its 
existence depends on the condition of irrotationality of 
the flow. For transonic flows over sufficiently thin 
bodies (small perturbation assumption!) it is usually 
possible to work with the velocity potential since any 
vorticity introduced by occurring shocks is of at least 
third order in the perturbation velocities. It follows 
that the velocity can be expressed as 



ii 

9X 



U = = <i> 



(n-6) 



V = = 



d± 
9 Y 






The continuity equation (II-2) may be written 



1 Dp _ . 

p Dt $ xx $ yy’ 



(II-7 ) 



Recall that the velocity of sound is 



2 ,a?s 

a = W 

K s=constant 



(II-8) 
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and that 



D£ = d£ DP = DP 

Dt dP * Dt 2 Dt* 

a 



( II-9 ) 



Therefore the continuity equation (II-7) becomes 



$ + $ + 
xx yy 



P +$ P +P, ) = 0 . (11-10) 

2 x x yy t 



a p 



Integration of the Euler equations along a streamline gives 
the unsteady Bernoulli equation 



♦t + ir + 1 t m p(t) 



(ii-iD 



1 2 

where F(t) = for uniform parallel flow at infinity. 

Taking the partial of (11-11) with respect to t, one has 



Ip = JL_[$ + K_] . 

p t 9t L *t 2 J * 



(11-12) 



Inserting (II-3) and (11-12) for P' , P , and P, in (11-10) 

X y T3 

and utilizing (II-6) one obtains the following partial 
differential equation in 



2 2 
$ * 

(1 —_)$ + (1 ^r-)$ 

k & 2 ' xx a 2 ;? yy 



2$ $> 

2 xy 



2$ 



xt 



^ ^y l _ 

f V ” 72 " *tt " 0 



(H-13) 
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The usual way to simplify this equation is to neglect 
all product terms obtaining 



( 1-M ) 4» + $ 



xx 



yy 



M 

_ CO 
2 -. — $ , 
a xt 



2 tt 
a 

CO 



= 0 



( 11 - 1*0 



which is the well known unsteady linearized potential 
equation valid for the supersonic and subsonic regimes. 
However, for the case near M =1. such a "naive" lineari- 

9 CO 9 

zation cannot, in general, be carried out. Consider the 
following analysis for the steady case. 

Define the perturbation potential, 4> , such that 



4> = U X + <J>. 
00 T 



(II-15) 



For the steady case, the Bernoulli equation (11-11) is 




constant . 



( 11 - 16 ) 



Now since 




Cp T 




R T 



and 

a 2 = (§£) = YRT 

9p s 



Equation (11-16) becomes 
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constant 



(II-17) 





The constant is evaluated by noting that at the critical 
Mach number 



W* = a* 

and (11-17) becomes 

a 2 = a* 2 - W 2 (11-18) 

With (11-18) and (11-15) in (11-13), one has 

A (a 2 -U 2 ) + A a 2 = A C(y+l)U <b d> 2 

v xx v yy w v xx LV ' °o^ x 2 ^x 

+ A ] + 2A CU A + A A ] 

2 Y y y xy L oo^y T x y y 

+ 4> yy [(Y-l)U co 4> x +^4> x 2 + *+! 4> y 2 ! (H-19) 

As was mentioned previously, under the Prandtl assump- 
tion, one would neglect the product terms on the right-hand 
side (RHS). However, for M^ = 1, a more careful analysis 
is in order. To this end one introduces the Prandt 1-Glauert 
transformation: 
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4> = e4>(*,Y) 



( 11 - 20 ) 



where X = X 

Y = Y l/l-M 2 ' 
00 

u 

00 

M = — 

« a 



( 11 - 21 ) 

( 11 - 22 ) 



Here, one requires <j> to be of order unity, therefore e 
assumes the order of magnitude of <p. 

Equation (II-18) becomes 



T ( i m 2 ^ 2 x T (1-M 2 )a 2 = e 2 (y+1 )U 4) — <j>- 
£<b — C 1-M )a + e<p — 03 <» 1 <» Y xx Y x 

' (Yi T i m r 



XX 



yy 






2 ^ li co / ^xx^y 



+ 2 e 2 <p — (1-M 2 )U <p- + 2e 3 <f) — d>— 4>— (1-M 2 ) 
Y xy <» oo^y Y xy Y x Y y ~ 



+ e 2 ( 1-M 2 ) (y-l)U ?-?■ — + e 3 ^(l-M 2 ) 2 <f> — <j>- 
<*> 1 oo x yy ■ 2 oo ' Y yy Y y 



+ £ 3 Cy.1) . (1-M 2 )(J) — 4> — 
2 °o /Y yy^x 



( 11 - 23 ) 



Now the terms on the RHS may be neglected for suitably 

small e . To investigate the validity of such a measure, 

one may compare the largest term on the RHS with a term 

2 

from the LHS. This term is the first term of order e , the 

3 2 2 

remaining terms being of order e or e (1-M ). Consider 
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the ratio 



^(Y+l) <t>- 
a 2 (1-M 2 ) 

CO S CO ' 



( 11 - 2 * 1 ) 



It is clear that the nonlinear term becomes of higher 
order for the subsonic and supersonic cases but must be 
retained for transonic flows. This consideration can be 
generalized to unsteady transonic flows thus leading to 
the transonic small perturbation equation as the governing 
equation: 



M M 

[1-M-(Y+D ~ <f> ]4 + <f> -2 — <}> , <f> = 0 

00 U Y x JY xx Y yy a xt 2 Y tt 

00 u co 

00 

(11-25) 



A rigorous justification for this equation was first given 
by Lin, Reissner, and Tsien [Ref. 11] and by Landahl [Ref. 12]. 

A further simplification of this equation is possible if 
the transonic flow is sufficiently unsteady. Although this 
can be shown by mathematical order of magnitude considera- 
tions in these same references, it is instructive to consider 
the physical aspects of the transonic wave propagation 
problem. Figure 3 shows an acoustic source in transonic 
flow. It may be seen that the parts of the wave fronts 
moving upstream, referred to as receding waves, move away 
very slowly from the source. The downstream or advancing 
waves move with a velocity approximately equal to twice the 
speed of sound. 
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Figure 3. TRANSONIC WAVE PROPAGATION 



Assuming now a continuous spatial distribution of distur- 
bance sources, it is readily understood from Fig. 4 that the 
receding waves will tend to cancel provided the wavelength 
(given by 2ir|a-U|/w , where w is the circular frequency of 
the disturbance) is very short compared to body reference 
length, C. However, if the sources oscillate very slowly. 





The Resulting Wave 
Figure 4. WAVE' CANCELLATION 
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pressure waves of the same sign will have sufficient time to 
interact and possibly to build up to finite amplitudes. 
Obviously, these receding waves may create large phase lags 
between motion and pressure on the airfoil causing severe 
aeroelastic instability problems. 

We can draw the conclusion, therefore, that sufficiently 
unsteady transonic flows can be analyzed by linearized 
theory, Eqn. (11-1*1), whereas only mildly unsteady transonic 
flows must be investigated using Eqn. (11-25). The condition 
for linearization is generally stated as [Ref. 12] 



K >> 



1 1 ' - M 



LOCAL 



where |l - flgoCAL I the largest deviation of the Mach 
number from unity in the flow and K is the reduced frequency 
(K = wC/Uj. 



B. DEVELOPMENT OF THE BOUNDARY CONDITIONS 

To completely describe the problem, proper boundary and 
initial conditions must be specified. The most important 
ones are: 

1. The flow on the surface of the body must be tangent 
at any instant of time. 

2. Sommerfeld’s radiation condition must be satisfied 
at any instant of time, i.e., waves propagate away from 
sources of disturbance. 

Further, certain edge and wake conditions must be 
satisfied, e.g., smooth flow at trailing edges must be 
established (Kutta condition). 
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The general flow tangency condition can be stated as 
follows. If the equation of the surface of a body which is 
moving in a time dependent manner is given by 

F(X, Y, t ) = 0, (11-26) 

then the flov; tangency condition requires that the total 
derivative of the function F is zero, i.e., 

~=F. +UF + V F =0. (11-27) 

Dt t x y 



Consider an airfoil that is placed in a Cartesian coordinate 
system in such a way that if the airfoil is restricted to 
small amplitude motions out of the Y = 0 plane, then its 
upper and lower surface can be described by 



F u = Y - Y u (X,t) = 0 
F l = Y - Y L (X,t) = 0 



( 11 - 28 ) 



Therefore, from the boundary condition one has 



V = Y + U Y 

u t U x 


at Y = Y u (X,t) 

(11-29) 


V = Y 4- U Y l 


at Y = Y L (X,t) 


t X 





One now introduces the small perturbation assumption, 

<i> = U M X + 4> (11-15) 



where the disturbance velocity components are 



such that 



(II-31) 

(II-32) 



U « + u 

V = v 

and u,v << U 

9 CO 

In (11-29) the term, u Y , can then be neglected compared 

u x 

to the much larger term, U C0 Y U . The boundary conditions 

x 

then simplify to 

v = Y u + U ot Y u for Y = Y u (X,t) 

(II-33) 

and v = Y t + U ro Y T for Y = Y T (X,t). 

L t L x L . 

* 

For small amplitude oscillations, (11-33) can be further 
simplified by expressing the flow tangency condition in the 
plane Y = 0. Taylor series expansion gives 

v(X,Y u ,t) = v(X,0 + ,t) + Y u ^^ .. x°- . xL) - + 

and v(X,Y L ,t) = v(X,0",t) + Y L a ^ X ?° -^1 + 

Using the same argument as above, the product terms can 
be neglected so that the flow tangency condition can be 
stated as 

v = Y + U m Y for Y = 0 + 

u, 00 u„ 

t x (11-34) 

v = Y, + U m Y t for Y = 0" 

L t L x 

For a thin foil performing a time dependent motion about 
the Y = 0 plane, the airfoil thickness terms can be dis- 
regarded by similar arguments. The above conditions become 
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V 



for Y = 0 



(II-35) 



= ah + it 111 
at u » ax 



where h(X,t) describes the time dependent chord line of the 
blade . 

C. THE DERIVATION OF THE UNSTEADY EQUATIONS 
Define the following nondimensional terms. 



X = 



Y = 



t = 



X 

c 

Y 

C 

t U 



(II-36) 



<f> = 



v = 



U C 

CO 



V 

u„ 



h = ~ 



where C is the chord of the foil. Equations (11-25) 
and (11-35) may now be written in nondimensional form: 

[1 - M 2 - (y+DM 2 ^-]^ — • 2M 2 ?-t- - M 2 $rr- = 0 

l oo ' 1 c ° Y x JY xx Y yy °° Y xt <*> Y tt 



(11-37) 



and 



V = h F + h- 



( 11 - 38 ) 



Since the main interest of this thesis centers on the 
determination of the aeroelastic stability of a foil, only 
small perturbations from the steady state need be studied. 
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If one restricts oneself to such small amplitude perturba- 
tions, it is possible to write the perturbation potential as 
the sum of a steady state and a time dependent potential. 
Also, in the case of a harmonic oscillation, one may describe 
the time dependent quantities as the product of a complex 
amplitude function and a harmonic time factor. Therefore 



one defines the potentials <J>, and such that 

?(X,Y,t) = ^(X,?) + iKX,Y,K)e iKF (11-39) 

and the foil chord function and vertical velocity are 

h(X,t) = h 1 (X) + h(X,K)e lKt (11-40) 

v(X,t) = v 1 (X) + C(X,K)e 1K ^ (II-41) 

where K = Jp (11-42) 

oo 

It may be shown that (11-37) becomes 

+ 4^— = 0 C II— 1 * 3 ) 



which is the steady state equation, and the unsteady equation 
is 

_m2 (y+ 1 )9_ (4^-ifj-) + 

3 x 

-2iKM*ip- + ^ = 0 . (11-44) 

The boundary equation (11-38) becomes 

V 1 ( X) = h 1 — (X) (11-45) 

and v(X,K) = h-(X,K) + iKh(X,K) (11-46) 

A 
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On the foil surface equation (11-43) must satisfy the 
boundary equation (11-45); (11-44) must satisfy (11-46). 

For convenience, the overbar and hat villi now be dropped 
and the steady and unsteady equations with boundary equa- 
tions become 



Cl-Mf - M ~ (Y+1)< ^lx ]<i) lxx + hyy = 0 

subject to 

v-l = h lx , t at Y=0 (11-48) 

and 

[1-M 2 ]<J> - M 2 (y+l)|~(4> 1 ip ) + M 2 K 2 i4 

00 r xx 00 1 9x lx y x 00 r 

- 2iKM 2 iJ; x + ip yy = 0 (11-49) 

subject to 

v = h + iKh , at Y=0 . (11-50) 

A 

Note that (11-47) is a nonlinear equation in Equation 

(11-49) is linear in but the existence of the variable 
coefficient <|>^ x considerably complicates its solution. 

D. A BRIEF REVIEW OF POSSIBLE APPROACHES TO 
UNSTEADY TRANSONIC FLOW ANALYSIS 

The following is a presentation of various approaches to 
the solution of equation (11-49). 

1 . Method of Teipel and Hosokawa 

An approximation of equation (11-49) was proposed 
by Teipel [Ref. 9] and Hosokawa [Ref. 13] for oscillating 
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airfoils in sonic flow (M TO = 1). For steady sonic flow, 
Oswatitch [Ref. 8] introduced a linearized approximation by 
setting 

(y+l)<J> lxx = A = constant, X > 0 , (11-51) 

allowing (11-43) to be written in the easily solved form 



<J> 



iyy 



x <j> 



lx 



o . 



(H-52) 



Teipel and Hosokawa generalized this procedure to unsteady 
sonic flow by approximating equation (11-49) with (neglecting 
the term in to ) 

a X 

-X ip x + 4> yy - 2iKi|> x + K^ip - 0 . (11-53) 

Defining Q 

Q = X + 2iK , (11-54) 

one has from (11-46) 



t yy - Q i4 x + K 2 ip = 0 . (11-55) 

This equation allows the study of sonic flow past a 
single oscillating airfoil in unbounded flow over the entire 
frequency range. It contains the steady solution (Oswatitch 
approximation) as a special case. For large K (Refer to 
page 23) the term containing 4> lxx becomes negligibly small; 
therefore, retention of X at high reduced frequencies 
becomes unnecessary. Thus, when one takes X = 0, the 
Teipel-Hosokawa theory also contains the linearized theory. 
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2 . Local Linearization Technique for 
Steady Transonic Flow 



It was shown that the Oswatltch approximation of the 
steady transonic small perturbation equation is 



♦lyy - X +lx = 0 > 



X = constant > 0 



(II-52) 



The solution to this equation for sonic flow past a single 
airfoil is 



•q(x,Y) = - -L, / x 



V(S) 

vAtT -i / x-s 



XY* 

i|(X-S)‘ 



ds . 



(II-56) 



Of particular interest is the velocity on the profile which 
is given by 



* lx (X,0) = - -L; If / X 
XX /vA QX -1 



V(s) 

/x-s 



ds 



(H-57) 



Recalling the definition of A, (11-51), Spreiter argued that 
the original variable, 'f’lxx 5 could be reintroduced into 
equation (11-57) thus linearizing the equation locally. 

This leads to an ordinary nonlinear differential equation in 
which can be solved analytically. The agreement with 
experiment is remarkably good [Ref. 10]. Its application 
to the choked wind tunnel wall interference problem was 
given by Sandeman [Ref. 4]. It is expected that an extension 
of this technique to the unsteady case may prove fruitful. 

3. WKBJ - Method 

Rubbert [Ref. I 1 !] applied the .WKBJ method to this 
problem. Differentiating the axisymmetric counterpart of 
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equation (11-^9) for = 1, multiplying by and applying 

the Hankel transform, the equation can be reduced to the 
form 

Y xx + = 0 (H-58) 

where N is a large parameter. Unfortunately no inversion 
of the resulting transform integrals has been found. 

4 . Method of Geometrical Acoustics 

A third approach has recently been proposed to obtain 
the thickness effects for oscillating wings in transonic 
flow [Ref. 15]. This procedure begins with a given non- 
uniform steady transonic flow field which has been deter- 
mined either analytically or experimentally (In the latter 
case an analytical expression must be obtained which matches 
the measured pressure coefficient.). Superimposing upon 
this flow field small unsteady perturbations, the path and 
transmission time of these acoustic perturbations are 
determined by the methods of geometrical acoustics. Again, 
only preliminary calculations have been carried out to date. 

5 . Unsteady Transonic Flow with Shock Waves 
Typically, transonic flows around an airfoil exhibit 

locally supersonic regions which are terminated by shock 
waves. Often, the interactions of small unsteady perturba- 
tions with these shock waves may be responsible for aero- 
dynamically induced instabilities. Eckhaus [Ref. 1 6 ] made 
the problem mathematically tractable by assuming a normal 
shock extending from the surface to infinity with uniform 
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steady supersonic flow in front and uniform steady subsonic 
flow behind the shock. Performing a small perturbation 
analysis, i.e., superimposing small time dependent perturba- 
tions upon this idealised steady transonic flow, he was able 
to obtain equations for the interaction of the sound waves 
with the shock. However, only a few numerical results were 
obtained . 

For the following study, the Teipel-Hosokawa 
approach was used. 
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III. TRANSONIC UNSTEADY CASCADE FLOW PROBLEM 



A. PROBLEM FORMULATION FOR SONIC FLOW PAST 
OSCILLATING UNSTAGGERED CASCADES 




2b *- 

Figure 5. UNSTAGGERED CASCADE 



Consider an infinite cascade of thin blades arranged as 
shown in Figure 5. Blade chord is 2b and blade spacing is 
P. Subscript zero referes to the reference blade and 
subscript one refers to the adjacent blade in the positive 
ordinate (n) direction. The reference coordinate system is 
fixed at mid-chord with the abscissa (£) in the direction 
of flow. 

Two n coordinates are specified with n Q originating on 
the reference blade and rij, on the adjacent blade such that 

= n 0 + P • (III-D 

Now the following nondimensional variables are defined: 

X = | (III-2 ) 
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(III-3) 



Y m = 



‘rn 

b 



171 = 0,1 



R 



P 

b 



(III- 1 !) 



As the governing equation, the Teipel approximation is 
used , 

^ yy - Q i|» x + K 2 = 0 . (11-55) 

This equation must be solved subject to the following 
boundary equations: 

1. Zero disturbance upstream of the cascade. 

2. Plow tangency on the blades. 

As shown in Section II-B, this latter requirement can be 
expressed mathematically as follows. Along the reference 
blade equation (11-50) may be written 

yx.V°> = , o = r* 1Kh 0 ’ (illI ' 5) 

and along the adjacent blade 

3 h, 

iP y (X,Y o =R) = Vl = + iKh 1 . (III-6) 

To solve this boundary value problem, the potential 
function if) is written as the summation 

iMX,Y o ) = $°(X,Y q ) + ^ 0 (X,Y o ) + *°(X,Y 1 ) 

+ $ 1 (X,Y 1 ) (III-7 ) 

where <J>° and are the isolated blade solutions for the 
o l 

reference and adjacent blades and <2> o and are interference 
potentials. 
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Such an approach has previously been proposed for the 
subsonic and supersonic cases by D. N. Gorelov [Refs. 3 and 
17]. Also refer to Samoylovich [Ref. 18]. 

Inserting (III-7) into (III-5) and (III- 6 ) and realizing 



C (X ’ 0) ' v m • ” ,= 0 > 1 > 



one obtains for the reference blade. 



*oy (X ’V 0) + *?y (X ’ Y l = - R > + ^iy (x > Y i = -R) = 0 



iy 



iy 



(ni-8) 



and for the adjacent blade. 



*oy (X ’V R) + $ oy (X,Y o = R) + V X *V 0) = 0 * 



iy 



(III-9 ) 



The isolated blade solution is already known from 
Teipel [Ref. 9] as 0 



2 

x exp[>Z(x-s) - £ 5 ^] 



v/x-S' 



- dS 



(III-10) 



and 



*?(X,Y 1 ) = 



VttQ 1 



x EXP[!^(X-S) - 9 5 ^] 

/ V,(S) — as 

-1 x y/x-S 

(III-ll) 



The interference potentials must have the same form but 

contain unknown source distributions u and u, : thus, 

o 1 5 ? ’ 

2 Y d 

x EXP[^(X-S) - | 

(X, Y ) = — f u„ (S ) 2 ■■■■■-- dS 

° 0 /iQ 1 -1 ° v^X-S^ 

( III-12 ) 
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and 



V X ’ Y 1 } = f u 

x x /irQ^ -1 



x' EXP[^r(X-S) 

1 / u-. (S) Q 



Q Li 

T x-s J 



v/x-s‘ 



dS. (in-13) 



These source distributions are approximated by the power 
series as follows: 



L-l 

u = y v .. x 

m n £ 0 mN 



N 



( 111 - 1*0 



and 



L2 



v = £ 0 .. 
m n £ 0 mN 



r N 



(HI-15) 



Inserting these approximations, (111-1*0 and (III-15), 
into equations (III-8) and (III-9) results in 



L-l 
N 



p on xN + 1 /f <w^N>/^^ Exp L lx - s) -?iL ]d; l 



= o (in-16) 



for the reference blade and 




V N+ 



R 

7 




,N 



(X-S) 






EXP[ 




(X-S)-| 




= 0 (III-17) 



for the adjacent blade. 
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Figure 6. BOUNDARY CONDITIONS FOR TT PHASE LAG 



The coefficients in (III-15) may be evaluated for the 
case of simple oscillation, no translation, and tt phase lag 
between blades as follows. One arbitrarily sets the angle 
of maximum amplitude for the reference blade at a Q = -1 
[Refer to Fig.. 6]. Then the angle of the adjacent blade is 
cx^ = 1. For small a, equation (III-5) may be written 

i|> (X, Y q = 0) = -1 - iKX (III-18) 

and (III-6) , 

^ y (X,Y 1 = 0) = 1 + iKX . (HI-19) 
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Hence, the coefficients are 



and 



G 

oo 



-1 



ol 



lo 



-iK 

1 




iK . 



( III-20) 



Similarly, for in phase oscillation, the coefficients are 



and 



6 = 0 . = -1 
oo lo 



6 , = 9,, = -iK . 
ol 11 



(I 11-21) 



Satisfying boundary equations (III-16) and (III-17) at 
L points along the blades generates a set of 2L equations 
in 2L unknown interference coefficients. The coefficients 
in hand, ^(X,Y o ) in (III-7) becomes a known function. 

A computer program for the solution of this problem in 
the tt phase shift case is developed in Appendix A. Some 
sample results are presented in Table II for the steady case 
and Table YII for the unsteady case. 



B. EVALUATION OP THE INTEGRALS OCCURRING 
IN EQUATIONS (III-16) and (III-17) 

The integrals occurring in (III-16) and (III-17) are of 
the form 

= exp[ 4 (x - s > - v c^st ] dS • (III - 22) 

By making the substitution, U = X-S, one obtains 

h = / X+1 - —y / ) - N - EXP[^ U - | 2i] dU (111-23) 

0 u 2 y 
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Recalling Q = X + 2iK, 1^ is complex and must be resolved 
into real and imaginary parts for evaluation. 

The exponential term in (III-23) may be written 



EXP[-q U 



$ ~] = EXP[ 



xk 2 u 
X 2 + 4k 2 



XR - 

Tuf- 



• EXP 



{ 1[ - 



2K 3 U 
X 2 + 4 K 2 



KR 

2U 



4 



( 111 — 21 } ) 



One now defines the following terms: 



A = 



r 



X 2 + ilK 2 



B = 



XR 2 

T~ 




Equation (111-2*0 may be rewritten, 



EXP[ib=r U 

y 



Q H 2 -, 

■ ir i r J 



EXP [ XAU - |]EXP{i[-2KAU - |] } 



(HI-25) 



Separating (III-25) into its real and imaginary parts 
obtains the following two fundamental integrals from 
(III-23) : 

|] C0S[2KAU + §] dU 
§•] SIN[2KAU + |] dU 

Where X N = J NR " 1J NI 



Xt-X/y TTNJN 

J NR = / J t EXP [ XAU - 

WK 0 u 



x+1, 



vN 



and I NI = / — j>" EXPCXAU - 

1MJ- Q U * 



, one 



(III-26) 



(III-27) 

(III-28) 
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Expanding (X-U) N in the finite binomial series, 



(x-u) M = l 



N+1 (-lP' 1 N! X N - J+1 U J - J 



J=1 



( J-l ) ! (N-J+l ) ! 



(HI-29) 



equations (III-26) and (III-27) become 



N+1 



k J-l 



.N-J+l X+l 



T _ V (~1) U ~‘ L N! x 1 " u ' rx f'LJ-f f , B-, 

X NR '(J-l)"! (N-J+l)! J 1 U EXP[XAU - y] 



• C0S[2KAU + |j]J dU 

(HI-30) 



and 



NI 



N+1 

1 

J=1 



(- 1 ) 



J-l 



tj-tt 



N! X N J+1 
(N-J+l) ! 



X+l 

/ 



u J -i" 



EXPCXAU - |] 



• SIN[2KAU + |]| dU 

(III-31) 

The fundamental integrals in (III-30) and (III-31) do 
not lend themselves to analytical integration. Numerical 
integration techniques may be used but care must be exer- 
cised in the vicinity of the lower limit. Here the integrand 
is singular when J < 5/2. The problem is further aggravated 
by the rapidly oscillatory character of the functions 
SIN[ AAU - jj] and COS[AAU - |] . 

Three specific cases are posed for the treatment of 
these integrals. 

1. The Steady Case 

Setting K = 0 reduces (III-30) and (III-31) to the 



following : 



and 



I NR 



N+l 



I 

J=1 



(-1) J_1 N! X N " J+1 
( J-l ) ! CN-J+1)! 




(III-32) 



I 



NI 



0 



The numerical integration of this case is straight- 
forward. The exponential term forces the integrand to zero 
near the origin permitting one to neglect that portion of 
the integral very near U = 0. Equation (III-32) may also be 
evaluated by expressing it in terms of the error function. 
(Note: complete agreement between the two integration 

techniques has been obtained thus verifying the numerical 
integration procedure for the steady case.) 

2 . The Oscillatory Case with 1=0 

In the unsteady case for X = 0, equations (III-30) 
and (III-31) contain the basic integrals 

X+l , r v 

/ U * C0S[2KAU + £■] dU (III-33) 

0 u 

and j +1 U J "^ SIN[2KAU + |] dU. ( III-3 4) 

0 



To avoid the previously mentioned difficulties near 
the lower limit, the integration must be accomplished in 
two steps , i . e . , 



J 1 U J ~* g^[2KAU + |]dU = / a U J * 2j I S j[2KAU + |]dU + 

+ / V"* g°®[2KAU + |]dU. 

a 

( HI-35) 

For judicious selection of (a), the second integral 
in the RHS may be numerically integrated in a straightforward 
manner. The integral on the limits 0 to (a), which will 
henceforth be referred to as the correction integral, may be 
approximated for sufficiently small (a) as follows. 

Recall the following relations: 

COS [ 2KAU + |] = C0S[2KAU]C0S[|] - SIN[2KAU]SIN[|] (111-36) 

and 

SIN[2KAU + |] = SIN[2KAU]C0S[|] + COS [2KAU]SIN[|] (III-37) 

and the truncated series expansions 

2 

C0S[2KAU] = 1 - ( 2 - K A u . ) . - (111-38) 

and 

SIN[2KAU] = 2KAU - i^KAU)_ . (III-39) 

Equations (111-38) and (III-39) are valid for 2KAU suffi- 
ciently close to zero. One may now write 



*\2 



✓ T ^ F 

f U * C 0 S[ 2 KAU + £]dU = 

0 U 



3 T * F 3 *r £ 

/ U * COS[=f]dU - ( 2KA ) / U * SIN[~-]dU 
0 u 0 U 

_ !2KA)_ j yJ-j C0S [E-j dU 
d 0 u 

+ . ( .. 2 . K . y . 3 . / a u J+ i S IN[|]dU (III-40) 

D o U 



and 



a T £ r* 

/ U * SIN[ 2KAU + £]dU = 



3 5" 3 ? 

/ U J " 2 SIN[|]dU + (2KA) / U J “*COS[|]dU 
0 u 0 U 

- /“ o J - * siH[£]an 

d 0 u 

(2KA) 3 ( a U J+S COS[|]dU . (III-41) 

5"- 6 U 



In theory integrals of the form 



a 



/ 

0 



t m/2 COS 
SIN 



[|]dU 



may be expressed in terms of the Fresnel integrals. Hence, 
integration of (III-40) and (III-41) can be accomplished in 
"closed" form. (In actuality there is a computer accuracy 
limit on this technique as described in Appendix A-IV. ) 



43 



Verification of this technique for J = 1 was accomplished 
by checking it against a known analytical function as 
follows : 

s?h C2kau + § ]du + [-^r s?n [2kau + l ]du - 
00 

/ 75 - s?N [2KAU + ! ]dU an-02) 

where the correction integral was evaluated using the Fresnel 
integral technique, the second integral was evaluated using 
numerical integration techniques, and the integral on the 
RHS was evaluated using the following relation: 

f & ^ ^ ] . (in-., 3) 

For the case K = 0.01 and a = 0.01, Table IA presents 
the comparison of the value of (III-43) with that obtained 
from (111-42) for various R's. (The numerical integration 
was truncated at an upper limit of 500 and the convergence 
criterion EPS for subroutine ZSUM was set at 1.0E-4.) 

Table IB displays the computed integral using (III-35) 
with an upper limit of X+l = 2.0 for various values of (a). 
(The convergence limit placed on ZSUM was set at 1.0E-5.) 

Table variables are listed below. 

DK = K 
DL = X 
RH0 = R 
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AL = a (Correction integral upper limit.) 



SCOR « 


= Value of correction integral (111-41) for J = 1. 


CCOR = 


= Value of correction integral (III-40) for J = 1. 


SUMI = 


= Value of numerical integration term from (a) to 
specified upper limit for the SIM integral. 


SUMR = 


= Value of numerical integration term from (a) to 
specified upper limit for the COS integral. 


SINT = 


- Value of (III-35) SIN integral for specified upper 
limit and J = 1. 


CINT * 


= Value of (III-35) COS integral for specified upper 
limit and J = 1. 


SSIN = 


= Value of (111-43) SIN integral. 


CCOS = 


= Value of (111-43) COS integral. 



3 . The General Case, K ^ 0 and 1^0 

In the apparent absence of a suitable correction 
integral for this case, one must evaluate the full equations 
(III-26) and (III-27) by numerical methods. The difficulty 
close to the lower limit is alleviated to some extent by the 
exponential term. The exponential function drives the 
amplitude of the integrand to zero more rapidly than a 
negative power of U can amplify the highly oscillatory 
behavior of the transcendental functions. For this reason, 
as in the case of K = 0, one may justify neglecting that 



portion 


of the integrals very near zero. 

Although no results are presented, the computer 


routine 


developed in Appendix A is designed to carry out 



such a case. 
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TABLE 



IA VERIFICATION OF CORRECTION INTEGRAL 



R 


SCOR 


SUMI 


SINT 


SSIN 


0.5 


34.62 


1. 016 


35. 63 


35.57 


1.0 


14.45 


3.45 8 


17.91 


17.85 


1.5 


4.968 


7. 032 


12. 00 


11.94 


2. 0 


-1.114 


10.16 


9.045 


8.987 


R 


CCOR 


SUMR 


CINT 


CCOS 


0.5 


15.48 


19.80 


35. 28 


35.32 


1.0 


-1.782 


19.34 , 


17.56 


17.60 


1.5 


-5.797 


17 .45 


11.65 


11.69 


2.0 


-4.490 


13.19 


8. 696 


8.736 



COMPARISON OF VALUES COMPUTEO USING 
EQUATION (I 11-43) WITH THOSE 06TAI NED 
USING (III —42 ) FOR AN UPPER LIMIT OF 500. 
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TABLE IB VERIFICATION OF CORRECTION INTEGRAL 



DK= 0.500 1 


DL = 


O 

• 

o 




RH0= 0.500 


AL= 0.010 


S COR= 1 . 57005 D 


00 




CC0R = 


1.671 14D-0 1 




S UM I = 3.96851D 


00 




SUMR = 


3 . 240 1 3 D 00 




S I NT = 5.53856D 


00 




C I NT = 


3.40724D 00 




DK= 0.500 1 


DL = 


0. 0 




RH0= 0.500 


AL= 0.050 


SCOR= 1.66484D 


00 




CC0R = 


-2.65483D 00 




SUMI = 3. 87 372 D 


00 




SUMR = 


6 . 0 6 2 0 7 D 00 




S I NT= 5. 53856D 


00 




C INT = 


3.40724D 00 




DK= 0.500 


DL = 


0.0 




RH0= 0.500 


AL= 0.100 


SCOR= 3. 74310D 


00 




CCOR = 


-1. 1280 ID 00 




SUMI = 1.79545D 


00 




SUMR = 


4.53526D 00 




S I NT= 5. 53856D 


00 




C I NT = 


3 . 4 0 7 2 4 D 00 




DK = 0.500 


DL = 


O 

• 

O 




RH0= 0.500 


AL= 0.200 


SCOR= 4. 62026D 


00 




CC0R = 


4. 962 7 0D— 01 




SUM I = 9.18297D 


-01 




SUMR = 


2.91097D 00 




SI NT= 5. 53856D 


00 




C INT = 


3.40724D 00 




DK= 0.500 


DL = 


0.0 




RH0= 0.500 


AL= 0.300 


SCOR= 4.87826D 


00 




CC0R = 


1 .27508 D 00 




SUMI = 6.60294D 


-01 




SUMR = 


2. 1321 6D 00 




S I NT = 5. 53856D 


00 




C INT = 


3.40724D 00 
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C. CALCULATION OP THE POTENTIAL FUNCTION AND 

THE PRESSURE COEFFICIENT 

The potential function is obtained from equation (III-7), 

*(X,Y 0 ) = *°(X,Y 0 ) + V X > Y o } + $ 1 (X » Y 1 ) 4 V X > Y i) 

(HI-7 ) 

The interference potentials and may now be 
computed using the previously determined coefficients (P 
m = 0,1; N = 0,...,L-1) in equations (III-12) and (III-13). 
Equation (III-7) may be rewritten in terms of the coeffi- 
cients as follows: 



♦ <X *V ■ J JJk ei 4' 5 - s >1 ife** 8 



, X -N f 2 o Y? 

+ •=; t®,/'™) / J^EXPC^tX-S)-^ ijxrsjldS 



1N I * 1 -l/xls 1 



(III-M) 



The resulting integrals may be evaluated using the same 
techniques utilized in the coefficient derivation. 

The pressure coefficient is obtained from the. relation 



Cp = -2 (4> x + iKi{0 . (111-45) 

By placing the partial with respect to X of equation (III-44) 
into (111-45), one may obtain the pressure distribution at 
any point in the flow field as well as on the blade surfaces. 
Further, the force and moment coefficients are easily 
obtained by integration over the blade chord. 



D. EVALUATION OP THE STEADY STATE 
POTENTIAL ALONG THE- BLADE 



The computation of the steady state potential along the 
upper surface of the reference blade follows from equation 
(III-44) evaluated at K = 0, Y q = 0 } and Y 1 = -R. The 
development which follows is for the ir phase shift case 
where [Prom Eqns. (III-20).] 



6 = -1 
oo 



6 lo = 1 



(111-46) 



Equation (III-44) becomes 

X 



wx,0) 



L-l 



X EXP[ - £■ ^-] 

/ _L_ dS + / ds 

-l /x-s‘ -i 



+ -f=, I 

/ttX n=o 



v 



X S N 



I VF5 EXP[ - T 5&s ]ds 



IN 



-1 /xls 



,N 



- V 



oN 



-l /x^s 1 



dS 



(111-47) 



The computer program is documented in Appendix B-I. The 
resulting potential values are compared with those obtained 
from the Sandeman wall interference solution of Section III-E 
in Tables III and IV. 



E. LAPLACE TRANSFORM SOLUTION FOR THE SONIC WIND 
TUNNEL WALL INTERFERENCE PROBLEM 

A parallel development to the collocation technique is 
found in the Laplace transform solution of the sonic wind 
tunnel wall interference problem [Ref. 4]. 



Virtual Blade 



R = 2d 




Figure 7 shows the arrangement of a blade facing a solid 
wall. Due to symmetry, the solution for the flow between 
the blade and the wall will be the same a^s that which would 
exist between the blade and the virtual blade beyond the 
wall. Hence it is the same solution expected of an unstag- 
gered cascade with it phase shift between the blades. 

The basic equation for oscillatory sonic flow from 

Section II is 2 

- Q4> y + = 0 . (11-49) 

j y A 



50 



The boundary-value problem can be formulated as follows: 

1. There is no disturbance upstream of the blade 
leading edge. 

2. Plow tangency conditions must be satisfied along 
the blade and along the solid wall. 

This latter condition may be stated mathematically. 





^y (x, o ) = 


at Y = 0 


(III-48 ) 


and 


<4 y (X,d) = 0 


at Y = ±d 


(111-49) 



Denoting the Laplace transformation of a function with 
respect to X by an overbar, i.e., 

CO 

f(p) =*£f(X) = / e“ pX f (X) dX , 

0 

one obtains for the Laplace transform of Eqn. (11-49) 

ip - X 2 i p = 0 (III-50) 

v yy r 

where X 2 = Qp - Y 2 (III-51) 

The boundary conditions transform into 

? = -v(p) Y = 0 (HI-52) 

and ip = ° Y = ±d (III-53) 

A solution to (III-50) satisfying (III-52) and (III-53) 
has been given by Miles [Ref. 19] > 

?(p,Y) = i v(p)CSCH(Xd)COSH[X(d-|Y| )]sgn(Y) ( III-5 1 ! ) 
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Restricting oneself _to the flow between the upper tunnel 
wall Y = d and the blade one has Y _> 0. Hence the absolute 
value sign may be dropped. 

To invert (111-5^0, a suitable expansion of the equation 
was given by Miles as 



?(p,Y) = j v (p ) 



oo 



I 

N=0 



(EXP[-X ( 2Nd+Y ) ] + EXP(-X[ (2N+2)d-Y] } ) 

(III-55) 



Recalling the Laplace transform of 



f(X) = X- 1/2 e- a/HX 



Rea > 0 



being 

00 Y 

f(p) = / e pA f(X) dX - /P , Re p > 0 , (111-56) 

0 . /p 1 



one has for 



f (p , Y) = j EXP[-X(2Nd+Y ) ] 



where 




the inversion 2 Q (2Nd+Y) 2 n 

EXPL-q X - -jr Y - 1 

f(X,Y) = 2 

/ttQX 1 



(III-57 ) 



Therefore, with the aid of the convolution theorem, one 
obtains the following solution: 



if>(X,Y) 




■ ; V -- S -T EXP[^(X-S)] • 

/xTs Q 

Q (2Nd+Y) 2 -, , rv?r Q [ ( 2N+2 )d+Y] 2 
f — x^s~ J + expl t X=S 



]j dS . 
( 111 - 58 ) 
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On the blade surface where Y = 0, this result simplifies to 



iKX,0) 



I 

/uQ' -1 



~ , EXP ['4 (X-S)] • 

/x^ y 



CO 




EXP[- 



2 2 
QN_d_ 

X-S 



] + EXP[- 



Q (2N+2) 2 d 2 
7 X-S 



3 



dS . 



J (III-59) 



For steady flov; where K = 0, one obtains the further 
simplification, 

f 

X 



<Kx,o) = - 



i 



V^rT 



/ 

-1 



v(S) 

/x-s' 



dS 



oo 



+ 2 I 

N=1 



f v(S) 



m 2.2 

EXP[-X — 33 ] dS 



(HI-60) 



This expression was previously given by Sandeman [Ref. 
in his analysis of choked wind tunnel flov;. Equation 
( 111 - 58 ) represents the generalization to oscillatory flow. 

Appendix B-II describes the translation of (III-60) into 
a computer program. The solution .of this function provides 
the steady state check of the collocation technique as is 
shown in Tables III and IV. 
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IV. PRESENTATION OF RESULTS 



Table II presents the interference coefficients for 
the steady state case obtained for the controlling param- 
eters indicated in the tables with the computer routine 
developed in Appendix A. 

Table III contains the value of the perturbation 
potential along the reference blade obtained from the 
routine developed in Appendix B-I and utilizing the inter- 
ference coefficients of Table II. These results may be 
compared with those in Table IV obtained from the steady 
state Laplace transform solution for the sonic wind tunnel 
wall interference problem of Section III-E and programmed 
in Appendix B-II. 

Table V demonstrates the erroneous correction integral 
obtained from Subroutine COR of Appendix A as the result 
of computer numerical inaccuracies. A comparison of the 
interference coefficients for two cases before and after 
the anomaly was detected is presented in Table VI. 

A sample of the unsteady interference coefficients for 
the case X = 0 is presented in Table VII. 

Table VIII compares a sample steady state output with 
the output of the general case X f 0, K / 0, for very 



small K. 



TABLE II INTERFERENCE FUNCTION COEFFICIENTS 
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TABLE II INTERFERENCE FUNCTION COEFFICIENTS 
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POTENTIAL SOLUTION FOR THE FOLLOWING PARAMETERS 
DL= 0.200 DK= C.O RHO= 1.000 

THE INTERFERENCE COEFFICIENTS 



DEL( 00) =-2. 045311 OE 00 
DEL(01)=-1.2504177E 00 
DEL {02 ) = 3. 0618846E--01 
DEL ( 03) =-l. 5100145E-G1 
DEL ( 04 ) = 1.0638958E-01 
DEL (05) =-6. 7890167 E-02 
DEL ( C6) =- 2 . 05 72 83 0E-02 
D EL { 07 ) = 3. 3797 889 E— 03 
DEL ( 08 ) = 1 . 3067 639E-0 1 
DEL ( 09) =-9. 7604275E — 02 



DEL (10)= 2.0453110E 00 
DEL ( 1 1) = 1. 2504177E 00 
DEL ( 12)=-3. 061884 6E -01 
DEL ( 13) - 1 . 5 100 145E-01 
DEL ( 14) =-l. 063895 8E- 01 
DEL ( 1 5 ) = 6. 7890 1 67E- 02 
DEL ( 16) = 2.05 72830 E-02 
DEL ( 17) =-3. 37 97 889E-03 
DEL ( 18) =-1.306 763 9E-01 
DEL ( 1 9) = 9.7604275 E-02 



VALUE OF 
STATION 


THE POTENTIAL ALONG THE 
P0TENTI AL 


BLADE 


X = -0 . 90 


POT = 1 . 1780910E 


00 


X 

II 

1 

o 

• 

CD 

O 


POT = 2.1 757832 E 


00 


X 

II 

1 

o 

• 

o 


P0 T= 3. 1754112E 


00 


X=-0 . 60 


POT = 4 . 1755638 E 


00 


X=-0. 50 


P0T= 5. 1756020E 


00 


X 

II 

1 

o 

• 

O 


POT = 6 . 1755772E 


00 


X=— 0. 30 


POT = 7 . 1755638 E 


00 


X=— 0.20 


P0T= 8 . 1 75 5724E 


00 


O 
<— 1 

• 

0 

1 

II 

X 


POT = 9.1755791E 


00 


o 

0 

« 

o 

1 

11 

X 


POT - 1 . 01 755 87 E 


01 


X 

II 

o 

• 

o 


P0T= 1.117 5572E 


01 


X= 0.20 


P0T= 1 . 217 5570 E 


01 


X= 0.30 


P0T= 1.31 75570E 


01 


X 

II 

O 

• 

O 


P0T= 1 . 41 75561 E 


01 


X 

II 

o 

• 

U1 

o 


P0T= 1 .5175563 E 


01 


X= 0.60 


P0 T= 1.6175568E 


01 


X= 0.70 


POT = 1 .7 17 5552 E 


01 


X= 0.80 


P0T= 1. 8175476E 


01 


O 

O' 

• 

o 

II 

X 


POT = 1.9175308E 


01 


X 

It 

H-* 

• 

o 

o 


POT = 2 .0176529 E 


01 
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POTENTIAL SOLUTION FOR THE FOLLOWING PARAMETERS 
DL= 0.200 DK= 0.0 RHO= 2.000 

THE INTERFERENCE COEFFICIENTS 



DEL (00) =-8. 0229187E-01 
DEL ( 0 1 ) =- 6. I 007 094E-01 
D E L ( 0 2 ) = 1 . 4 25 3 49 7 E- 01 
DE L( 03) =-6.6 098273 E-02 
DEL ( 04 ) = 2.7 520720E-02 
DEL ( 05 ) =— 1 . 2785759E-02 
DE L ( 06) = 5. 68 7 223 7 E-02 
DEL ( 07)=- 5. 71 69359E-02 
DEL (08) =-4.0 132940 E-02 
DE L ( 09) = 5. 057951 8E-02 



DEL ( 10)= 8 . 02291 8 7E- 0 1 
DE L ( 1 1 ) = 6 . 1007 094 E-0 1 
DEL ( 12) =-l. 42 53 A 9 7E — 01 
DEL (13)= 6. 6098273E — 02 
DEL(14)= -2. 7520720 E-02 
DEL ( 15)= 1. 2785759E-02 
DEL ( 16) =-5.68 72237 E-02 
DE L ( 17) = 5. 71 69 3 59 E-02 
DEL (18)= 4. 0132940E-02 
DEL (19) =-5.05 79 5 18 E-02 



VALUE OF THE POTENTIAL ALONG THE BLADE 



STATION 


POTENTIAL 




o 

• 

0 

1 

II 

X 


P0T= 8 . 1 58 14 14E- 


•01 


o 

00 

• 

0 

1 

II 

X 


POT = 1 .3204842E 


00 


X 

II 

1 

o 

• 

-sj 

o 


P0T= 1 . 8202248E 


00 


o 

o 

• 

0 

1 

11 

X 


POT = 2.3200750E 


00 


X=-0. 50 


POT = 2 . 820 0855 E 


00 


X 

It 

1 

o 

• 

o 


P0T= 3.320 1075E 


00 


X =-0 . 30 


POT = 3 .8 20 1075E 


00 


X=— 0. 20 


P0T= 4 . 320 0874 E 


00 


X=-0.10 


P0T= 4.82007 50E 


00 


O 

o 

• 

0 

1 

11 

X 


P0T= 5 .3200960 E 


00 


o 
< — 1 

• 

o 

II 

X 


P0T= 5.8200798E 


00 


X= 0.20 


P0T= 6 . 3200665E 


00 


X= 0.30 


P0T= 6 . 8200712 E 


00 


X 

II 

o 

• 

O 


P0T= 7.3200731E 


00 


X = 0.50 


POT = 7 .8200769E 


00 


X= 0.60 


P0T= 8. 3200636E 


00 


X= 0.70 


POT = 8 . 8200598E 


00 


X= 0.80 


POT = 9 .320 1199 E 


00 


X= 0.90 


P0T= 9. 8202019E 


00 


X= 1.00 


POT = 1.03 1 9364E 


01 
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POTENTIAL SOLUTION FOR THE FOLLOWING PARAMETERS 
DL = 0* 200 DK= 0.0 RHO = 5.000 

THE INTERFERENCE COEFFICIENTS 



DEL ( 00)=-l . 153845 8E-01 
DEL (01 ) =-l . 8925 136 E-01 
D EL ( 021= 7. 1354136E — 03 
DE L ( 03 ) = 2. 2944096E-02 
DEL (04) =-2.5 75 1 479 E -02 
DEL ( 05 ) = 2. 8709877E— 02 
DEL(06)=-6. 2098 347 E— 02 
DE L ( 07) = 4. 1337617E — 02 
DEL ( 08 ) = 5. 2630220E-02 
DEL ( 09 ) =— 5 . 1 547259 E— 02 



DE L ( 10)= 1. 1538458E-01 
DEL (11)= 1 . 8925 136E-01 
DEL ( 12) =— 7. 1354136E-03 
DEL ( 13) =- 2 . 2944096E-02 
DEL ( 1 4 ) = 2.575 1 479 E-02 
DEL ( 1 5 ) =-2 « 8709877 E-02 
DEL (16)= 6. 2098 347E — 02 
DEL ( 17) =-4. 1337 617 E-02 
DEL ( 18) =-5. 2630220E-02 
DEL ( 19) = 5 . 1547 259E- 02 



VALUE OF THE POTENTIAL ALONG THE BLADE 



ST AT ION 


POTENTIAL 




X=- 0. 90 


P0T= 7. 9726553 E- 


•01 


X = - 0 . 8 0 


P0T= 1 . 1 28 2368E 


00 


X 

It 

1 

o 

« 

o 


POT = 1. 38 5 5448 E 


00 


X 

II 

1 

o 

• 

o 

o 


P0T= 1.61151 7 0E 


00 


X=-0 . 50 


POT = 1 . 8 232851 E 


00 


X 

II 

1 

o 

e 

o 


POT = 2 . 02 861 4 0 E 


00 


X=-0. 30 


P0T= 2.231 0429E 


00 


O 

CNJ 

• 

0 

1 

II 

X 


P0T= 2 .432 1480E 


00 


o 

rH 

• 

0 

1 

II 

X 


P0T= 2 . 632 6466 E 


00 


o 

o 

• 

0 

1 

II 

X 


POT = 2 . 8 328 724E 


00 


X= 0.10 


POT = 3 .0329742 E 


00 


X= 0.20 


P0T= 3. 2330236E 


00 


X 

II 

o 

• 

U) 

o 


POT = 3 . 4330444E 


00 


X 

II 

o 

• 

o 


P0T= 3 .633 0500 E 


00 


X 

II 

o 

• 

o 


P0T= 3 . 8330574E 


00 


X = 0.60 


POT = 4 .0 330648 E 


00 


X= 0.70 


P0T= 4.233 0627 E 


00 


o 

CO 

« 

o 

II 

X 


POT = 4.4330502E 


00 


X= 0.90 


POT = 4 .6330080 E 


00 


X= 1.00 


P0T= 4 . 8337231 E 


00 
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POTENTIAL SOLUTION TOR THE FOLLOWING PARAMETERS 
DL= 0.200 DK= 0.0 RHO= 10.000 

THE INTERFERENCE COEFFICIENTS 



DEL ( 00) =- 1. 5652659E-03 
DEL (01)=-8 . 5002407E-03 
DELI 02 ) =-l .4888730E-02 
DELI 03 ) =— 5. 3 168200E — 03 
DELI 04 ) = 7.090T287E-03 
DELI 05) =--6. 0203 695E— 04 
DEL ( 06)=- 4. 930614 5E- 03 
DEL 107) = 3. 8378909E— 03 
DE L I 08) = 1. 6905249E-03 
DELI 09) = - 2. 18 47 31 9E- 03 



DE L 1 1 0) = 1. 5652659E-03 
DELI 11) = 8. 5002407E — 03 
DEL I 12 ) = l . 438 873 OE- 02 
DE L I 1 3 ) = 5.3168200E-03 
DELI 14) =- 7. 0907207E-03 
DEL 1 1 5 ) = 6 „ 028 3695 E-04 
DELI 16)= 4. 9836145E— 03 
DEL I 17) =-3. 8378909E — 03 
DEL 1 18) =-1.69 85249 E-03 
DELI 19)= 2. 18473 1 9E-03 



VALUE OF THE POTENTIAL ALONG THE BLADE 
STATION POTENTIAL 



x= 


-0. 


,90 


POT = 


7 


.9786545 E- 


-01 


x=- 


-0. 


, 80 


POT = 


1 


. 1 28 37 03 E 


00 


X =■ 


-0. 


.70 


POT = 


1 


. 38 1 9685 E 


00 


x= 


-0. 


, 60 


POT = 


1 


.5957613 E 


00 


x= 


-0. 


, 50 


POT= 


1 


* 7841234E 


00 


X = 


-0. 


.40 


POT = 


1 


.9544506E 


00 


x= 


-0. 


30 


POT = 


2 


.11 11937 E 


00 


x= 


-0. 


.20 


POT= 


2 


. 2573242E 


00 


x= 


-0. 


,10 


POT = 


2 


.3949909 E 


00 


x= 


-0. 


. 00 


POT = 


2 


. 5258188E 


00 


x= 


0 . 


.10 


POT = 


2 


.6510839E 


00 


x= 


0. 


, 20 


POT = 


2 


.7718019E 


00 


x= 


0. 


.30 


POT= 


9 

4 — 


. 8888006E 


00 


X = 


0. 


.40 


POT= 


3 


.00275 14E 


00 


x= 


0. 


. 50 


POT = 


3 


. 1 1 420 06 E 


00 


x= 


0, 


.60 


POT= 


3 


. 2235994E 


00 


x= 


0. 


.70 


POT = 


3 


.3 3 1 31 22 E 


00 


x= 


c. 


.80 


POT = 


3 


. 4376440E 


00 


X = 


0 , 


.90 


POT = 


3 


. 542 8400E 


00 


x= 


1. 


, 00 


POT = 


3 


.647 13 76 E 


00 
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TABLE IV VALUE OF THE POTENTIAL ALONG THE BLADE — WALL INTERFERENCE TECHNIQUE 
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TABLE IV VALUE OF THE POTENTIAL ALONG THE BLADE — WALL INTERFERENCE TECHNIQUE 



o 

o 

o 

« 

o 



II 

o 

X 

o'. 



oo 



Ql 






t-H 


O 


O 


O 


O 


O 


O 


O 


0 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


UJ 






O 

1 


O 


O 


O 


O 


O 


O 


O 


0 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


UJ 






1 

UJ 


LU 


UJ 


U) 


UJ 


UJ 


UJ 


UJ 


UJ 


UJ 


IU 


U! 


LU 


LU 


LU 


UJ 


IU 


UJ 


UJ 


UJ 








LPi 


00 


00 


r- 


CO 


vO 


O 


CO 


0 


CM 


O 


r-A 


O 


in 


O 


O 


CM 


in 


in 


r~A 


< 






00 


co 


0 > 


I s - 


f — 1 




CM 


CO 


0 


CO 


O 


00 


CO 


r- 


CM 


O 


CO 


v 0 


CO 


-H 


CL 






00 


00 


t-A 


tn 




> 4 - 


r~A 


r- 


in 


in 


rH 


r~C 


CO 


CM 


M* 


CO 


r-H 


I s - 


CM 


r- 


< 




_j 


I s - 


(M 


CO 


o> 


co 


LO 


r-A 


in 


o> 


CM 


m 


h- 


00 


O 




CM 


CO 


CO 


4 ' 


M* 


Ci. 


0 


< 


0 > 


r— I 


CO 


LT\ 


r- 


O 


rH 


CM 


CO 


in 


0 


I s - 


00 


O 


r-H 


CM 


CO 


> 4 - 


in 


0 


O 


# 

0 


H- 


I s - 


r-H 


1 — I 


H 


r-A 


r-H 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


co 


CO 


CO 


CO 


CO 


C0 


co 
















































* — < 




UJ 










































:s 


11 


H 










































0 


x 


O 










































_j 


0 


CL 


11 


II 


II 


II 


II 


II 


II 


II 


11 


II 


II 


II 


II 


11 


II 


II 


II 


II 


11 


II 


-j 






h- 


1 - 


h- 


h" 


H 


h- 


H- 


h- 


h- 


l~ 


H* 


\- 


h~ 


H- 


H 


H* 


h- 


H- 


h- 


b- 


0 






0 


O 


O 


O 


O 


O 


0 


0 


0 


O 


O 


CJ 


O 


O 


O 


O 


O 


O 


O 


0 


U- 






Cl 


CL 


CL 


CL 


CL 


CL 


Cl 


CL 


Cl 


CL 


CL 


O. 


CL 


CL 


Cl 


Cl 


CL 


CL 


CL 


CL 



UJ 

X 



cC 

o 

LL O 

o 

z CM 

o • 

►H O 

H 

ID 

-J II 
O -J 
oo O 



2 


O 


O 


0 


0 


0 


0 


0 


O 


0 


0 


0 


0. 


0 


0 


O 


0 


0 


0 


0 


0 


O 


O' 


CO 


I s - 


vO 


in 




CO 


CM 


» — 1 


0 


H 


CM 


CO 


M* 


in 


0 


I s - 


co 


o> 


a 


H* 


O 

1 


0 

1 


0 

1 


0 

1 


0 

1 


0 


0 


O 

a 


0 


0 

1 


0 


O 


O 


O 


O 


0 


0 


0 


0 


1 — i 


H* 


1 

II 


1 

II 


1 

II 


1 

II 


I 

II 


l 

n 


1 

II 


1 

11 


I 

II 


1 

II 


II 


II 


II 


II 


II 


11 


II 


II 


II 


11 


CD 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 



64 



(\J 

o 

o 

o 

o 

in 

CM 











DC 


00 




































• 


< 


1 — 




































• 




» — t 






























< 






c£ 




21 


3 


rH 


CO 


O 


CO 


0 


CM 


in 


O 


vO 


in 


Mb 


in 


Mb 


> 






< 




r-i 


< 


O 


0 


O 


0 


rH 


rH 


rH 


rH 


rH 


T-~i 


rH 


rH 


rH 












3 


DC 


1 


1 


1 


1 


1 


1 


1 


J 


f 


! 


1 


1 


1 


w 






b- 






O 


O 


0 


O 


0 


O 


O 


O 


O 


O 


O 


O 


O 


0 












Z! 


LU 


in 


CM 


<r 


b 


0> 


in 


in 


CO 


Mb 


CO 


O' 


rH 


a> 


PQ 


O 




5l 




O 


b- 


00 


rH 




CO 


0 


rH 


in 


0 


rH 


O' 


CM 


b 


M" 


< 


O 










Z! 


00 


rH 


in 


O 


CO 


in 


Mb 


00 


vO 


in 


00 




Mb 


E-< 


in 










*— » 


rH 


O 


Cb 


CO 


in 


st* 


in 


rH 


CM 


b* 


vO 


b 


CO 




« 








X 
































0 




a: 


O 


X 


00 


rH 


rH 


0 


vb 


CM 


rH 


b 


rH 


CO 


CM 


(M 


CO 


CM 








LU 


O 


LU 


O 


















1 






1 


1 




11 




a. 




w(M 


0 






























O 




Cl 


O 


00 O 
































X 




3 


in 


O | 
































o: 






r- 1 


OO 






























00 






CC 


00 


O 






























c£ 






0 


00 


rrun 






























UJ 






LL 


Mb 


CJ(M 






























h~ 








• 


• 






























LLJ 






UJ 


0 


ZX) 




































3 




1 — \ 






























< 






3 




00 II 






























C£ 






< 


11 


— UJ 






























< 






> 


0 


in 






























Ou 








0 


































0 






0 


CMO 


3 


O 


CO 


in 


b 


0 


CM 


Mb 


O 


O 


in 


Mb 


in 


Mb 


o 


0 




< 




1 LL 


< 


O 


0 


0 


O 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


rH 


2: 


in 




CC 




3 


DC 




1 


1 


1 


1 


1 


1 


I 


1 


| 


1 


1 


| 


M 


• 




O 




wO 


O 


O 


0 


0 


a 


O 


O 


O 


O 


O 


O 


O 


O 


O 


15 


O 




LU 




* 2 : 


LU 


00 


in 


0 


M* 


Mb 


b 


O 


rH 


CO 


CO 


00 


M" 


O' 


O 






b* 




Vr < 


b- 


0 


b 


0 


O' 


O 


co 


CO 


in 


0 


rH 


CO 


ro 


b 


3 


II 




z: 




X 


z: 


CM 


b 


in 


0 


rH 


0 


00 


b 


00 


CO 


ro 


in 


CM 


3 


X 


O 


>-< 




v-ro 


► — < 


rH 


<b 


O 


CM 


O 


00 


C 0 


O' 


O' 


in 


HO 


rH 




O 


O 


UJ 






30 






























LL 




GO 


3 




< 1 






in 


CM 


.-h 


vO 


CM 


rH 


rH 


rH 


Mb 


rH 


M 3 


rH 






3 


LU 




cca 


►— < 






















1 


1 




LU 








0 


00 


00 




























X 




z: 


00 


0 


LUO 






























b- 




0 


LU 




HO 


































*— < 


(X 


O 


z: • 






























cL 




h- 


LL 


CM 
































O 




0 




CM 
































LL 




UJ 


OO 


00 


a 11 


































a: 


O 


00 


1 — — 1 






























OO 




a: 


O 


CO 


< 






























2 




O 




• 


z 






























O 


O 


O 


O 


0 


00 






























►H 


• 




Z 




Hh 






























b- 


0 


•> 


< 




b- 






























3 




0 




II 


30 






























3 


11 


11 




00 


3 oC 






CM 


CO 


<b 


in 


O 


b 


00 


O' 


0 


H 


CM 


CO 


O 


3 


3 


► — < 


00 


□ UJ 


3 




















3 


rH 


rH 


rH 


00 


O 


O 


00 


00 


00 M 































65 



SOLUTIONS FOR THE FOLLOWING PARAMETERS TABLE VB 
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INTERFERENCE FUNCTION COEFFICIENTS 
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TABLE VIA INTERFERENCE COEFFICIENTS BEFORE AND AFTER ERROR DETECTION 
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TABLE VIB INTERFERENCE COEFFICIENTS BEFORE AND AFTER ERROR DETECTION 



TABLE VII INTERFERENCE FUNCTION COEFFICIENTS 



o 








pw 






PW 


p-% 


p~* 


pw 


o 






PW 


p> 


p-> 


p-% 




PW 


PS 


o 


o 


rH 


r—A 


r-A 


H 


(XI 


CM 


CM 


r-A 


*H 


o 


pH 


r-A 


r—A 


CM 


rH 


rH 


rH 


rH 


rH 


A 


a 


o 


o 


o 


o 


o 


O 


O 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 




1 


1 


1 


1 


1 


1 


1 


1 


1 


r-A 


1 


1 


1 


1 


1 


1 


1 


1 


1 


o 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


o 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LLI 


• 


<M 


O' 


r— 1 


O' 


I s - 


A 


r- 




O 


(XI 


• 


O' 


A 


A 


CM 


A 


h- 


I s - 


CO 


O 


o 


r- 


rH 


r- 


A 


o 


A 


r — 1 


h 


A 


O 


o 


O' 


O 




O' 


CM 


O 


CM 


o 


A 




o 


vt 


in 


r- 


CM 


CM 


o 


O' 


r-A 






A 


A 


A 


r- 


O 


(XJ 


vO 


rH 


r- 


II 


O' 


r- 


<1- 


00 


CM 


I s - 


r—A 


s0 


I s - 


»H 


II 


O' 


rH 


CM 


v0 


-t 


(XJ 


o 


(XJ 


>4* 












































< 


1— 1 

1 


I s - 

i 


>4- 


r—A 

1 


r—J 


O 


A 

i 


r-A 

I 


r—A 


r—A 

1 


< 


00 

1 


A 

■ 


(XJ 


O' 

i 


rH 


rH 

1 


A 


A 


> 4 * 




1 


i 

•. 


► 


1 

•» 




1 

•* 


i 


1 


*+ 


1 

#- 




1 


1 

c> 


•» 


i 


*» 


1 

t» 


•» 


► 


*» 




rH 


r— t 


(XI 


(XI 


CM 


CM 


CM 


CM 


CM 


CM 




CM 


rH 


CM 


CM 


CM 


r—A 


rH 


rH 


CM 




o 


0 

1 


o 


o 


O 

1 


O 

i 


O 

i 


O 

I 


O 

i 


O 

l 




0 

1 


O 


O 


O 

■ 


O 


0 

1 


0 

1 


0 

1 


O 




LLJ 


1 

LU 


1 

UJ 


1 

LU 


1 

LU 


1 

UJ 


i 

LU 


1 

LU 


1 

LU 


i 

UJ 




1 

LU 


1 

LU 


1 

LU 


1 

LU 


1 

LU 


1 

LU 


1 

LU 


1 

LU 


1 

LU 




r—A 


r- 


r- 


(XI 


h- 


h- 


*4* 


O' 


00 


CM 




A 


CO 


rH 


O' 


O 


o 


vO 


vO 


O' 




o 


CM 


o 


O' 


O 


rH 


CM 


A 


A 


A 




O 


O 


O' 


v 4 * 


A 


A 


r- 


rH 


A 




o 


in 


vO 


00 


A 


00 


O' 


O 


O 


O' 




A 


CO 


> 4 * 


(XI 


CO 


A 


> 4 - 


> 4 * 


O 






a 


00 


O' 


00 


A 


r- 




vO 


r- 




H 


CM 


CO 


h- 


r- 


CM 




vO 


A 


o 


• 


• 


• 


t 


• 


• 


• 


• 


« 


• 


o 


• 


• 


t 


• 


• 


• 


• 


• 


t 


o 


vO 


O' 


vO 


C\J 


A 


A 


A 


vO 


O' 


O' 


o 


O' 


^ 4 * 


A 


rH 


A 


rH 


rH 


A 


O' 


in 










1 






1 


1 




o 








1 


1 






1 


1 


• 


w 


WP 


wp 


w 


wp 


WP 


wp 


w* 


WP 


WP 


* 


WP 


WP 


WP 


WP 


w* 


WP 


WP 


WP 


•wo 


o 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


r-H 


II 


II 


II 


II 


II 


II 


II 


II 


II 




P-Si 


p— * 


PW 


P'S 


P-* 




p«* 


p-^» 


*—» 


p-« 




pw 


*—* 


PW 


r-* 


p-w 


pw 


p^ 


p> 




II 


o 


rH 


(XI 


a 


* 4 - 


A 


o 


I s - 


00 


O' 


II 


o 


r-H 


CM 


A 


vt 


A 


o 


h- 


CO 


a 


rH 


rH 


rH 


» — < 


rH 


rH 1 


r-A 


r—A 


rH 


H 


o 


r—A 


rH 


rH 


rH 


rH 


rH 


r-A 


rH 


rH 


X 




wp 


w 


wp 


wp 


WP 


■wp 


w* 


WP 




X 


WP 


WP 


WP 


w 


WP 


WP 


WP 


WP 


WP 


a: 


-J 


—1 


.J 


_J 


_J 


—1 




_J 






OC 


_J 


-J 


_J 


—1 


_J 


—1 


-J 


-J 


-J 




LU 


LU 


LU 


LU 


LU 


UJ 


LU 


LU 


LU 


LU 




LU 


LU 


LU 


LU 


LU 


UJ 


LU 


IU 


LU 




O 


o 


O 


O 


o 


o 


O 


O 


O 


a 




O 


a 


O 


Q 


O 


o 


O 


o 


O 





r—m 


p« 


p-^ 


— S 


«PW 


pp> 


P-W 


PW 




PW 




PW 


—■S 


PW 


p-> 


-P-* 


PW 


p"s 


PW 


PW 




o 


rH 


rH 


rH 


rH 


CM 


CM 


CM 


r-A 


rH 




rH 


rH 


rH 


CM 


rH 


rH 


r—A 


rH 


rH 




o 


0 

1 


0 

1 


o 

1 


0 

1 


0 

1 


O 

I 


O 

i 


o 

1 


0 

1 




0 

1 


o 

1 


O 

1 


O 

t 


o 

1 


O 


o 

1 


0 

1 


0 

1 




LU 


1 

LU 


1 

UJ 


1 

UJ 


1 

LU 


1 

LU 


i 

LU 


1 

UJ 


1 

LU 


1 

LU 




1 

LU 


1 

UJ 


1 

LU 


1 

LU 


1 

UJ 


1 

LU 


1 

LU 


1 

LU 


1 

LU 


O 


CM 


O' 


rH 


O' 


h- 


A 


h- 


't 


o 


CM 


o 


O' 


A 


A 


CM 


A 


I s - 


r- 


00 


O 


O 


r- 


rH 


h- 


A 


O 


A 


rH 


rH 


A 


O 


o 


O' 


O 


> 4 - 


O' 


(XJ 


O 


CM 


O 


A 


A 


o 


st* 


A 


r- 


CM 


CM 


O 


O' 


tH 




A 


A 


A 


A 


h- 


o 


CM 


vO 


rH 


I s - 


« 


O' 


h- 


> 4 * 


00 


CM 


h- 


rH 


o 


r- 


rH 


• 


O' 


r—A 


CM 


O 


> 4 " 


(XJ 


O 


(XJ 


> 4 - 


O 






















O 


« 


• 


• 


• 


t 


« 


t 


t 


• 


|| 


rH 




Nt 

1 


rH 


rH 


vO 


A 


rH 


r-A 


rH 


|| 


00 


A 


CM 

1 


O' 


r—A 

| 


rH 


A 


A 

1 




II 






l 




1 












II 






1 




1 






i 


1 


o 


rH 


rH 


(XJ 


(XJ 


CM 


CM 


(XJ 


(XJ 


(XJ 


CM 


o 


CM 


rH 


CM 


(XJ 


CM 


rH 




rH 


CM 




0 

1 


0 

1 


o 


0 

1 


O 

i 


O 

1 


o 

1 


o 

1 


O 

i 


O 




O 


O 

i 


O 

i 


O 

■ 


O 

1 


O 


0 

1 


0 

1 


O 

■ 




1 

LU 


1 

LU 


LU 


1 

UJ 


1 

LU 


1 

LU 


1 

LU 


1 

IU 


i 

IU 


1 

LU 




1 

LU 


1 

LU 


1 

LU 


1 

UJ 


1 

LU 


LU 


i 

LU 


1 

LLI 


1 

UJ 




»H 


h- 


h» 


(XJ 


r- 


h- 


4 “ 


O' 


00 


CM 




A 


CO 


rH 


O' 


O 


vO 


vO 


s 0 


O' 




o 


CM 


O 


O' 


o 


rH 


(XJ 


A 


A 


A 




O 


o 


O' 


Mf 


A 


A 


I s - 


rH 


A 




o 


A 


vO 


00 


A 


00 


O' 


O 


vQ 


O' 




A 


00 


XT 


(XJ 


00 


A 


>4* 


■4- 


O 




*4- 


A 


00 


O' 


00 


A 


r- 


*4* 


sO 


n- 




rH 


CM 


CO 


r- 


r- 


CM 




vO 


CM 




• 


• 


• 


• 


• 


« 


• 


• 


• 


« 
























O 

1 


O' 


vO 

1 


(XJ 

1 


A 


A 

1 


A 

1 


vO 


O' 


O' 

i 




O' 

1 


*4* 

1 


A 

1 


rH 


A 


H 

I 


rH 

| 


A 


O' 




1 

WP 


w* 


1 


1 


WP 


i 


1 




WP 


1 




1 


1 


1 


WP 


WP 


1 

w* 




WP 


wr 




II 


II 


II 


II 


II 


II 


II 


II 


II 


ii 




II 


II 


II 


II 


II 


II 


II 


II 


II 


o 






p-^ 




PN 


pp<^ 








r> 


o 


rw 


P^S 






PW 


*p> 


— «* 






• 


o 


rH 


CM 


A 


vf 


A 


vO 


h- 


00 


O' 


• 


o 


rH 


CM 


A 


*4" 


A 


vO 


h- 


00 


o 


o 


O 


O 


o 


o 


O 


o 


O 


o 


o 


o 


o 


o 


O 


O 


O 


O 


o 


O 


o 




WP 


WP 


WP 


WP 


WP 


WP 


WP 


WP 


WP 


WP 








WP 


WP 


w* 




WP 


•w^ 


WP 


II 




-J 


— J 


-J 


-J 


—1 


_J 




— J 


-J 


II 


—I 


_J 




1 


_J 


-J 


_J 




-J 


_J 


UJ 


LU 


UJ 


LU 


UJ 


LU 


LU 


LU 


LU 


LU 


_J 


LU 


LU 


LU 


UJ 


LU 


LU 


LU 


LU 


LU 


o 


o 


O 


o 


O 


o 


O 


O 


o 


O 


O 


o 


O 


a 


O 


o 


o 


a 


O 


O 


O 



69 
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DEL(09)=(-3.7490E 01, 8.8355E 00) DEL(19)=( 3.7490E 01,-8. 8355E 00) 
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TABLE' VII INTERFERENCE FUNCTION COEFFICIENTS 
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TABLE VIII INTERFERENCE FUNCTION COEFFICIENTS 
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V. DISCUSSION OP RESULTS 



A computer program was designed to compute interference 
coefficients for the collocation technique developed in 
Section III-A. For the steady case, the resulting coeffi- 
cients in Table II were applied to the technique developed 
in Section III-C and III-D to provide the values of the 
blade surface potentials in Table III. For this same case, 
the sonic wind tunnel wall interference solution of Section 
III-E provided the data presented in Table IV which verified 
the collocation technique for that case. 

Although no "separate source" checks were made for the 
coefficients presented in Table VII (The case of simple 
oscillation about the mid-chord with no translation.), every 
effort was made to verify the procedure. Table VIII com- 
pares the proven steady state solution with a general small 
K solution. Unsteady integral approximation techniques were 
justified in Sections III-B-1, -2, and -3 as summarized 
below. 

Difficulties arising from the solution of the occurring 
integrals in the unsteady solution were not insurmountable 
but required careful handling. Whenever numerical integra- 
tion of functions involving ^ or (^) near X = 0 is 

attempted, accuracy is suspect. One must closely inspect 
arguments, limits, and results to ensure validity. For 
example, for the unsteady case and X ^ 0, with no apparent 
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means of computing the integral in the immediate vicinity 
of the origin, one must decide, against the shortcomings 
in computer capacity and programming methods, just what 
approximation must be applied and, once applied, how 
accurate it is. In this case, the question arose as to 
how closely one might approach the origin with numerical 
integration techniques before accuracy was degraded. At 
the same time the required computer time increased 
dramatically as the integrand became more ill-behaved. 

The most difficult question that arose in this investi- 
gation was the accuracy of the correction integral in the 
case x = 0. Once this approach w as found valid for J = 1 
through the calculations shown in Section III-B-2, there 
remained the need to improve the accuracy of the technique 
and ensure its validity for J > 1. In question were the 
number of terms needed in the truncated series shown in 
Equations (III-40) and (III-41), the validity of the inte- 
gral reduction formula used in subroutine RED, and the 
selection of the upper limit of the correction integral. 

Although four terms were used in the truncated series 
for the case J = 1 (Refer to Appendix A-IV), only three 
terms affected the significant figures of the results for 
the parameters used in this investigation. For J = 2, 3> 
and 4, the number of terms used was decreased by one for 
each computation. The primary motivation for this step was 
the inaccuracy of the integral reduction formula for large 
J and large R. 
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Referring to Tables ,VA and VB, it was determined that 
numerical inaccuracies began to occur for large J in the 
integral reduction routine. The effect of this inaccuracy 
was not important for small R but for the case of R = 5, 
for example, it overwhelmed the results. For integrals 
of the form 



^/ AL u J l COS (E/U )dU, 

as J becomes large and for AL << 1, the value of the integral 
must become small. However, in the case of R = 0.5, the 
magnitude of the integral values diminished with increasing 
J up to J = 9 but increased for J > 9* In the case of 
R = 5> the anomaly in the computation affected integrals 
for J’s as low as five. (In this latter case, it is obvious 
that the erroneous answer for J = 13 must affect the results.) 
The cause of this error was determined to be differences 
taken between large numbers in the reduction formula. The 
value obtained for the Fresnel integrals was only as accurate 
as the curve fit used in Subroutine CS. This error and the 
error generated in computing the SIN and COS functions 
occurring in the formula resulted in erroneous answers of 
a magnitude that was directly related to J and R. To 
correct the problem, the use of the integral reduction 
routine for J >_ 5 was eliminated. 

For J >_ 5, numerical techniques were resorted to. Here, 
for AL << 1, the amplitude of the integrand was small 
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permitting one to neglect that portion of the integral near 
zero. This technique seemed to produce reasonable results, 
at least the results did not have large magnitudes that 
would have affected the final answers. The resulting 
integral values were in some cases less than the convergence 
limits placed on the numerical integration subroutine, ZSUM, 
and were, therefore, assumed to contain error not greater 
than the said limits. Further, as these integral values 
became small, it should have been possible to neglect them 
altogether since the values of those integrals with small 
(J) predominated in the series. 

Comparative results for the case K = 0.5 and R = 0.5 
and 5-0 are presented before and after the error was 
detected in Tables VIA and VIB. 

The upper limit of the correction integral was selected 
such that the approximation ^g(E/U) = ^^(DU + E/U) is 
valid in the range 0 <_ U <_ AL. The point AL also needed to 
be within the capabilities of the numerical integration 
scheme on the limits AL <_ U <_ X + 1. Such a point was 
picked at AL = R/100 which ensures that (For U £ AL) E/U 
was at least l.OE^ larger than DU. Such an AL was also 
within the capabilities of the numerical technique over the 
range of parameters used in this investigation. Table IB 
demonstrates that the use of the truncated series in this 
approximation decreases its sensitivity to AL. 

Computer time could have been reduced with a larger AL 
if one were willing to accept reasonable inaccuracies as 
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the term DU became more significant. A smaller AL increased 
computer time and the accuracy of the numerical integration 
technique became more and more questionable as the integral 
became more ill-behaved. Computer time for a single run 
of the program (for X = 0 ) on the IBM 360/67 computer 
system was on the order of 15 minutes. 

It should be noted that in the absence of a good error 
measurement yardstick, statements concerning more or less 
error are indeed relative. It is difficult to assess 
where computer numerical error or approximation error' be- 
comes more important. It is for this reason that an inde- 
pendent check of these results against another source is 
necessary. The Laplace transform solution for the sonic 
unsteady wind tunnel wall interference problem would, in 
theory, be such a check; however, this technique will 
result in integrals of the same form as those used in the 
collocation solution and will require similar (or the same) 
solution techniques. It will be difficult, then, to obtain 
a completely independent check of these results. 

Once verification or modification of results presented 
in this investigation is in hand, application of Spreiter's 
local linearization technique alluded to in Section II 
should be applicable and lead to reasonable analytic 
results for transonic pressure and moment distribution 
prediction. 
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VI . SUMMARY 



Following the basic approach offered by Gorelov for 
subsonic and low supersonic flows, a collocation technique 
was applied to the purely transonic case superimposing the 
isolated foil solution provided by Teipel with unknown 
interference potentials. A computer routine was developed 
which would generate the coefficients for the unknown 
interference source distributions in the steady case. The 
interference coefficients in hand, the potential along the 
blade was determined through a second routine. These 
results were favorably compared with the potential solution 
derived from a routine based on the Laplace transform 
treatment of the sonic wind tunnel wall interference problem 
for the steady case. All three of these computer routines 
initially utilized series expansions which reduced the 
occurring integrals to functions of the well known error 
function. Confirmation of these results led to an extension 
of the collocation approach into the unsteady regime. 

First attempts to utilise series approximations for 
the unsteady case resulted in failure due to nonconvergence 
of the series. Faced with a lack of an apparently simple 
solution to the resulting integrals, numerical integration 
was resorted to. Here difficulties arose centering around 
the pathological behavior of the integrands near the origin 
as pointed out in the Discussion of Results. 
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As a basic numerical integration technique, the "pyramided" 
Simpson's rule method was incorporated. For the steady case 
where the integrals were easily computed, the results of 
this method were checked against the steady state case 
already in hand and found to agree perfectly. 

The greatest difficulty occurred in the unsteady case 
for X = 0 where one could not reasonably neglect the inte- 
grals in the interval near zero. To overcome this problem 
an analytical approximation (for the integral evaluated 
near zero) was introduced and checked against a known 
analytical solution shown in Section III. The technique 
involved the reduction of the integrals to functions of 
the Fresnel integrals. However, this procedure resulted in 
additional problems which were not immediately detected. 
Although the technique was theoretically sound, the computer 
numerical error pointed out in the Discussion of Results 
became apparent after study of the output. Once the cause 
was detected, the anomaly was remedied in a simple manner. 

The coefficients for the unsteady interference functions 
as presented here shall require a "separate source" check. 

The Laplace transform solution for the oscillatory transonic 
wind tunnel interference problem should provide such a 
check. If these techniques verify one another, a third 
independent check could be the Hamamoto linearized high 
frequency solution. 

Development of required computer routines for the 
potential and pressure coefficients from the collocation 
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and Laplace techniques should be facilitated by the use of 
subroutines derived for the unsteady coefficient routine 
presented here. 

It is felt that the collocation technique will prove a 
versatile tool in transonic flov; studies in that the free- 
dom which one has in specifying boundary conditions will 
permit treatment of problems which defy solution by other 
techniques. For example, the sonic wind tunnel wall inter- 
ference technique is only good in the it phase shift case 
while the collocation technique will provide results for 
any phase relationship. Further, it should be valid for 
more complex boundary conditions such as those posed by 
wind tunnels with porous walls. 
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APPENDIX A 



FORMULATION OF THE COLLOCATION COEFFICIENT ROUTINE 



A-I. MAIN PROGRAM 



The purpose of the main program is the solution of 
Equations (III- 16 ) and (III-17) for the unknoivn inter- 
ference function coefficients. 



Equation (III-16) must be satisfied for each station along 
the lower blade while simultaneously equation (III-17) must 
be satisfied at each station along the upper blade. For 
the purposes of this routine, ten stations have been selected 
(X = -0.9+0.2r, r = 0,1,2, ... ,9) • Computations near the 
leading edge are avoided due to the singularity there. 

For a given K, X, and R, the term 





/ — 

-1 (X-S) 



X 0 N 




(III-17) 




-1 (x-s) J/z 



X oN 

/ — 



; EXP[^-(X-S)-^ |ig]dS = SR(X,N)+iSM(X,N) 

( A-I-l ) 
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is a complex function of the parameters X and N. The 
complex value of this function at each (X,N) is provided 
by Subroutine WSUM and will be dealt with under Appendix A-II. 
WSUM returns the value of (A-I-l) as real and imaginary 
parts in arrays SR and SM respectively. 

The problem then is to evaluate the various known or 
computed values and insert them in proper form in a suitable 
matrix for simultaneous solution. 

Observe (III-16) and (III-17) become 



With 0j N and 0 q determined by boundary conditions at each 
of the ten stations, (A-I-2) and (A-I-3) form a series of 



Taking the set of coefficients, 0 MN , from (III-21) , 
(A-I-2) and (A-I-3) are rewritten for the in phase case as 




^ in )(SR(X,N) + iSM(X,N))J 



} - 0 

(A-I-2) 



and 




twenty equations in unknowns and N = 0,1,2,. ..,9. 

Note, for the unsteady case, and are complex. 




+ i(SM(X,0) + DK SR(X,1)) 



(A-I-4) 
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and 



1 {?: IN xN+ ^0l/ SR ^ X,N ^ + iSM(X}N ^j = (SR(XjO)-DK SM(X,1)) 
+ i(SM(X,0) + DK SR(X,1)) (A-I-5) 



where DK replaces K in the computer routine. 

For the ir phase shift case, the coefficients (III-20) 
apply and the lower blade boundary equation (A-I-2) 
becomes 



^2^P 0N X N +^ IN (SR(X,N) + iSM(X,N))J = -(SR(X,0)-DK SM(X,1)) 
- i(SM(X,0) + DK SR(X,1)). (A-I-6) 

The system of twenty equations may be written in matrix 

form 



X D = B 



(A-I-7 ) 



where X, D, and B are the matrices of coefficients, unknown 
/)' s, and right hand sides. X, D, and B are complex and may 
be written 



X = X x + X 2 i 

D = ^ + D 2 i ( A-I-8) 

and B = B^ + B 2 i. 
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Then in (A-I-7) 



(X x + X 2 i)(D 1 + D 2 i) = B 1 + B 2 i 
or (D-^-D^) + KDg^+D^) = ^ + Bgi. (A-I-9) 

It may be concluded that 



D 1 X 1- D 2 X 2 " B 1 



D 1 X 2 +D 2 X 1 = B 2 



or 



< 





( A-I-10 ) 



The twenty complex equations in twenty unknowns 
represented by (A-I-7) have been transformed into 40 real 
equations in ^0 real unknowns. One has only to solve 
for and D 2 and set 

P = D. + iD 2 . 

PROGRAM STEP ANALYSIS 

(Refer to Routine MAIN at the end of this Appendix. ) 

Double precision is used throughout the routine for 
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numerical accuracy. The use of the various arrays and 
common storage areas that are declared will be detailed 
below. 

Various arrays and fixed parameters are initialized 
in MAIN to minimize computational redundancy in later pro- 
gramming. The arrays COC(IO), COS (10), AS (10), and AC (10), 
zeroized in DO LOOP 100, are utilized in Subroutine V/SUM. 

DL = X, DK = K, and RHO = R, the controlling parameters, 
are initialized and PI = it is set. The value for 



is computed. For the case DK = 0, A is set to zero to 
avoid possible computer problems involving division by 
zero. Additional parameters are computed as follows: 



A 




( A-I-ll ) 



(DL) 2 + 4(DK) 2 



B 



DL(RHO) 2 

Ij 



( A-I-12 ) 



C 



DL(DK) 2 



(A-I-13) 



(DL) 2 + 4(DK) 2 



D 




(A-I-14) 



(DL) 2 + 4 ( DK ) 2 



DK(RHO) 

2 



2 



E 



(A-I-15) 



G = | TAN(^) 



(A-I-16) 
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Here G is set to 



G = tt/^ 



( A-I-l6a) 



for the case DL = 0 to avoid division by zero. 



AK = SIN(G) 



(A-I-17) 



BK = COS (G ) 



(A-I-18) 



Z = RHO [(DL) 2 + 4(DK) 2 ] 
2 v/p I 



(A-I-19) 



For future ease in changing the number of stations on 
the blades, the parameter L is used. Here L = 10 is set. 

The initial station, X^ = -0.9, and the interval between 
stations, XINT = 0.2, are set. 

The first write statement is included to indicate the 
controlling parameters DL, DK, and RHO. 

Finally, the matrix manipulation section begins. 
Referencing (A-I-10), array AAC^OjAO) is dimensioned to 
hold the coefficient matrix. BB(^O) holds the matrix of 
right hand sides before solution and the matrix of unknowns 
after solution. (A-I-10) reflects the diagonal duplication 
of the real and imaginary submatrices, X^ and X^ . Note 
that equations (A-I-^O and (A-I-5) form adjacent rows in 
the array and that there is duplication between them for a 
given X. The first two rows of X^ are as follows: 
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1 



SR(X ,0) SR(X ,1) . 



'X. 



SR(X 13 9) 



SR(X i5 0) SR(X 1 ,1) 



. . . SR(X 1 ,9) 



( A-I-20 ) 



corresponding to the unknowns 



^ROO ^ROl * • * ^R09 ^RIO ^Rll * ' ' ^R19 



where ^RMN + ^IMN* two rows of Xg are 



0 si(x 1 ,o) SI(X 1 ,1) ... si(x 1 ,9) 

0 ... 0 



si(x r o) si(x 1 ,i) ... si(x 1 ,9) o 



( A-I-21 ) 



corresponding to 

^100 ^101 



Ao 9 ?■ 



111 ^112 ^ 119 * 



Reference to (A-I-10) shows the placement of these rows. 
The first row of X^ fills AA(1,1) to AA(1,20) as well as 
AA( 21 , 21 ) to AA ( 21 , ^1 0 ) . The first row of X ? fills AA(1,21) 
to AAClj^O), and, with signs reversed, AA(21,1) to AA(21,20). 
For the ir phase shift case, the right hand side matrix is 
filled four rows at a time as follows: 
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BB(1 ) = 


-(SR(X 1 ,0) 


BB(2 ) = 


(SR(X 1 ,0) 


BB ( 21 ) = 


-(si(x 1 ,o) 


BB ( 22 ) = 


(si(x 1 ,o) 



and 



- DK SI(X ,1)) 

- DK SI(X 1 ,1)) 

DK SR(X 1> 1)) 

+ DK SR(X 1 , 1 ) ) 

( A-I-22 ) 



BB ( 3 ) 
BB(4) 
BB(23) 
BB( 24 ) 



(SR(X 25 0) - 

(sr(x 2 ,o) - 
(SI(X 2 ,0) + 
(SI(X 2 ,0) + 



DK SI(X 2 ,1)) 
DK SI(X 2 ,1)) 
DK SR(X 2 ,1)) 
DK SR(X 2 ,1)) 



and so forth. 



DO LOOP 1000 on J = 1,L' varies the station from X-j to 
X^q. It sets two rows in submatrices X^ and X 2 , or, 
equivalently, four rows in AA and four in BB. A is used 
to set the power of X in submatrix X^ and is initially set 
to A = X° = 1. JJ is a counter used to select every other 
row, i.e., JJ = 1,3. •••,19 as J = 1,2,. ..,10. 



With X^ set at X = -0.9, WSUM is called to provide the 
integral term in the boundary equations. WSUM computes 
the integral for X^ and for N = 0,1,2,. ..,9. These ten 
values are returned in common arrays SR and SM. SR holds 
the real parts and SM the imaginary as shown in functional 
notation in (A-I-l). The arrays are filled as follows: 
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SR( 1 ) = SR(X, 0 ) 



SM( 1 ) = SM ( X , 0 ) 



SR( 2 ) = SR(X, 1 ) 



SM ( 2 ) = SM ( X , 1 ) 



(A-I-23) 



SR(10) = SR(X, 9) 



SM (10) = SM ( X , 9 ) 



DO LOOP 900 on I = 1,10 controls the columns. The 
columns of submatrices X^ and X^ are filled in sequence for 
each station counted by DO LOOP 1000. In each step in 
LOOP 900, four positons in submatrix X^ are filled with the 
power of X. Pour positions in X^ are zeroized. The real 
part of the applicable integral is placed in four positions 
in X^ via the dummy variable SUMR. Similarly, the imaginary 
part is placed into Xg via SUMI. Care is taken to obtain 
the proper signs. Once the LHS is set for each station, 
four rows in the RHS matrix BB are set according to (A-I-22). 

The station is then augmented by the interval, XINT, and 
the procedure is repeated until all positions in AA and BB 
are filled. 

For verification, the matrix prior to solution is printed 
out. Matrix BB is destroyed in the solution process. 

The matrix is solved for the unknowns using library 
Subroutine DGELG. The inputs are the coefficient matrix, 

RHS matrix, coefficient matrix size (here set to ^0 by M =i IL), 

RHS matrix size, and a significance measure (here set to 1.0E-6). 
The solution is returned on BB and any error indication on IER. 
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BB(1) to BB(10) contain the real parts of the coefficients 
and BB(ll) to BB(20) contain the imaginary parts. The parts 
are combined in complex form and placed in complex array ZB. 
The form of ZB is 

ZB ( 1 ) = DELT (0,0) = 

ZB( 2 ) = DELT (0,1) = 



ZB(ll) = DELT ( 1 , 0 ) 



(A-I-24) 



ZB ( 20 ) = DELT (1,9) = 

Finally the desired coefficients are printed out on 
computer paper output and on cards. 
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A-II . SUBROUTINE WSUM 



The purpose of WSUM is to provide the control for the 
computation of the integral 



I 

e 



R 

2 



X S N K 2 O R 2 

— — rEXPC^-CX-S)-^ —q]dS, 

-i (x-s r 2 Q 1 X " S 



/ 



( A-II-1 ) 



As is demonstrated in Section III-B, (A-II-1) may be 
reduced to 



I 



e 





( A-II-2 ) 



where 



I NR 



N t 1 (-1) J “ 1 N!X N “ J+1 
j^U-1) ! (N-J+l) ! 



X+l 

/ 

0 



(' 



u J "/ 



EXP[ AAU-B/U ] 



(II-30) 



• COS [2KAU+E/U] f dU 



and 



I 



NI 



N+l 

E 

J=1 



( -1 ) J-1 N ! X N J+1 
(J-1) ! (N-J+l ) ! 



x +!r j_£ 

/ fu 2 EXP[ AAU-B/U ] 

0 1 

(II-31) 



• SIN [2KAU+E/U] J dU 



where the fixed parameters A 3 B 3 C 3 D 3 and E are set in MAIN 
[Refer to Eqns. (A-I-ll, -12, -13, -1*1, and -15).] and 
carried to WSUM in common space HOLD. 
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Given the integral values, the summation is a simple 
computer manipulation. The solution of these integrals 
involves three variations of the same basic numerical inte- 
gration technique. The numerical integration of (III-30) 
and (III-31) is accomplished by a "pyramided" Simpson's 
rule using end correction. This technique is presented in 
Appendix A-III where Subroutine ZSUM is described. ZSUM 
prepares the numerical solution of the basic integral for 
whatever values of J,K,A, and upper and lower limits, XO 
and XA, are specified. 

The difficulty encountered with the numerical integration 
of this basic integral occurs as the result of the lower 
limit, XA = 0. For J < 2.5, the amplitude of the Integrand 
may be unbounded near XA = 0. Further, except in the special 
case where K = 0, the SIN and COS functions in the integrand 
become unmanageable near zero where the controlling term 
in their arguments, E/U, becomes large. The resulting 
infinitesimal period cannot be treated by numerical tech- 
niques. For a reasonable A, one may neglect the troublesome 
portion of the integral near zero. This is permissable 

since the exponential function of (-B/U) forces the ampli- 
/ 

tude to zero in the interval where the integrand is ill- 
behaved. When A = 0, for J < 2.5, one may not neglect the 
integral near zero. Fortunately, in this case, one may 
approximate the ill-behaved portion of the integral by 
analytical methods. 
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For A = 0, the basip integrals in Eqns. (III-30) and 



(III-31) become 



X+l 



,J-f 



SUMR = / U * COS [DU + E/U]dU (III-33) 

0 



and 



SUMI = f 
0 



Y±1 T— ^ 

U *SIN[DU + E/U]dU (III-3 J 0 



Introducing the correction term of Section III-B, one may 
write ( III-33) and (III-S^) as 



X+l t _£ 

SUMR = COC + . U * COS [DU + E/U]dU (A-II-3) 

AL 



and 



X+l T r 

SUMI = COS + f U *SIN[DU + E/U]dU ( A-II-4 ) 
AL 



where the correction terms are 



AL j_r 

COC = / U * COS [DU + E/U]dU 

0 



(A-II-5) 



and 



al T_r 

COS = / U 2 SIN [DU + E/U]dU 

0 



(A-II-6) 



Integrals COC and COS are prepared by Subroutine COR, 
The technique as presented in Appendix A-V provides an 
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approximation which reduces the correction integral to the 
well known Fresnel integrals. 

The limit AL is selected such that (1) the numerical 
integration technique may be applied to the second term 
in (A-II-3) and (A-II-4), and (2) the correction approximation 
is valid over its interval. 

STEP ANALYSIS 

(Refer to the Subroutine at the end of this Appendix.) 

Subroutine WSUM is called to provide the block of L 
integrals at each station on the blades where the boundary 
conditions are enforced. (In this investigation L is set 
at ten.) The subroutine inputs are DK = K, DL = A, RHO = R, 
and X. 

In the first steps several fixed parameters are estab- 
lished for convenience. EPS = 1.0E-8 is the convergence 
parameter used in the numerical integration process in 
Subroutine ZSUM. LIM = 50 is a step limiter also used by 
ZSUM [Refer to A-III]. PI = tt is used in later calculations 
in WSUM. 

A test is made to determine which of the three cases 
is being dealt with. If K = 0, branching is executed to 
Step 1. If not, the next step is encountered. If K / 0 
and A = 0, the routine branches to step 2. The case K / 0 
and A / 0 follows at step 3. 

Step 1, K = 0 

The intention here is to use the numerical integration 
process over the entire interval. Since no SIN or COS 
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functions enter, one nay, approach closely to zero at the 
lower limit (XA is set to 1.0E-8). The upper limit is set 
to XO = X + 1 . Subroutine ZSUM is called to compute the 
integrals J = 1,2,..., 10. The integral in (11-30) is 
returned on SUMR and the negative of the integral in (11-31) 
is returned on SUMI. (Note that in this case SUMI = 0.) 

The values are stored in the vector arrays AC and AS. The 
routine then jumps to step 1000 where the summation 
indicated in Eqns. (11-30) and (11-31) takes place. 

Step 2, 1 = 0, K ? 0 

In this step the special case of X = 0 is computed. 

The Fresnel integral is obtained, and where applicable, the 
correction is computed and summed with that portion of the 
integral obtained by numerical methods. 

Observing that the value of X does not affect the 
correction integrals, they need be computed only once. 

Hence a test is made: 

IF(I.GT.l) GO TO 22. 

Step 22 bypasses the correction computation. 

The first steps in computing the correction integrals 
includes the preparation of inputs for Subroutine COR and 
Subroutine Red. The rationale for these functions is given 
in A-IV and A-V. They include the following: 



AL = 


RHO/100 


(A-II-7) 


AR = 


E/AL 


(A-II-8) 



105 



SS = SIN(E/AL) 


(A-II-9) 


CC = COS (E/AL) 


( A-II-10 ) 



The Fresnel integrals for upper limit AH are obtained 
from library Subroutine CS. The explanation for CS is 
given on call from the library and a copy is found at the 
end of this Appendix. The only argument required is AR. 
The following is returned. 



CCC 



SSS 



AR 




COS (T) 
/2ttT' 







SIN ( T) 
/2uT’ 



dT 



dT. 



(A-II-lla) 

(A-II-llb) 



In preparation for their use in Subroutines COR and RED, 
(A-II-lla,b) are modified as follows: 



CCC = [| - CCC] (A-II-12a) 

» 

and 

SSS = [~ - SSS], - (A-II-12b ) 

Refer to A-V for an explanation of these functions. 

The various inputs in hand. Subroutine COR is called to 
prepare the block of L = 10 correction integrals for Eqns . 
(A-II-5, -6). The arguments of COR are AL, D, E, CC , SS, 
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V 



CCC, and SSS. The valup of (A-II-5,-6) for each J is 
returned on CCOR and SCOR respectively and stored in arrays 
COC and COS. COC and COS are of vector size ten and are held 
in common with MAIM for permanency. Also note that (-SCOR) 
is stored in COR for compatibility with the numerical solu- 
tion of the remaining integral. The correction terms are 
printed for reference. 

Finally the remainder of the integrals are computed 
numerically in the same manner as was used in step 1 but 
with lower limit AL. The result is summed with its correction 
and stored in AC and AS. 

Step 3. K 0, 1 £ 0 

In this step, one relies on the exponential function to 
reduce the amplitude of the integrand near zero. That portion 
of the integral from zero to 1.0E-3 is neglected. The 
numerical solution is called as before and the results 
stored in AC and AS. 

Once the values of the integrals for J = 1,2,..., 10 are 
in hand, the indicated summations, Eqns. (111-30,-31), 
must be carried out. Step 1000 begins this operation. The 
summation is carried out in a straight forward manner for 
each of the L = 10 integrals. 1^ and 1^ are computed on 
SUMC and SUMS respectively. (Note that IFAC(M) is a function 
subroutine generating the factorial of M. Refer to A-VI.) 

The completed integrals are stored in arrays SR and SM and 
are printed for reference. 
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Recalling that equation (A-II-2) contained Q = 
a final complex manipulation must be made. Observe 

Q = \/l 2 + hK 2 [COS (to) + iSIN(w)] 

where 

to = ARCTAN ( 2K/X ) . 



Then 



= Z[BK + iAK] 

where 

, . R h 2 + UK 2 ] * 

2/tT 

AK = SIN(G) 



and 



BK = COS ( G ) 



for 

G = | ARCTAN ( 2K/A ) . 



Finally, (A-II-2) may be written 



X B = Z(I NR BK " I NI M) + 1Z(I NR AK + 



A+2iK, 

that 

(A-II-13) 



(A-II-14 ) 

(A-I-19) 

(A-I-17) 

(A-I-18) 

(A-I-16) 

(A-II-15) 
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Z, AKj and BK are computed in MAIN to eliminate 
repetition and are carried to VISUM in common. DO LOOP 1030 
computes the real and imaginary parts of (A-II-15) for the 
ten integrals and returns these values to MAIN in arrays 
SR and SM. 
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PERFORM FINAL IMAGINARY OPERATION ON INTEGRAL VALUES 
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A-III. SUBROUTINE ZSUM 



Subroutine ZSUM computes the value of the integral in 
the general Eqns. (III-30) and (III-31) by numerical 
methods . 



SUMR = / u * EXP[CU-B/U]COS[DU+E/U]dU (A-III-1) 

XA 



xo T_r 

SUMI = f V° * EXP[CU-B/U]SIN [DU+E/U ]dU (A-III-2) 
XA 



The technique utilized is Simpson's rule with end correction. 

Noting the behavior of the integrand, especially near 
XA = 0, it is not satisfactory to compute the integral using 
a fixed interval. Instead a "pyramided" scheme is used 
which automatically reduces the size of the subintervals 
until the scheme can compute the integral to within the 
limits of a satisfactory convergence criterion. 

The general form of Simpson's rule with end correction 
is 




SUM = Y5(7Y 1 +l6Y 2 +l4Y 3 +l6Y iJ + . . .+16Y N _ 1 +7Y N ) 




t I 



or 



N-l N-3 



SUM = 7 ( Y, +Y M )+l6 SY 2i +l4z Y 

15 1 N 1=1^ i=l 




(A-III-3) 
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where and are the value of the Integrand and the 

derivative of the integrand evaluated at station (i), and 
H is the length of the subinterval [Refer to Pig. 8]. 




Figure 8. INTEGRATION DIAGRAM 



Figure 9 demonstrates the "pyramid" algorithm. X 1 and 

X^ are the first and last station of the interval over which 

Y(X) is to be integrated, i.e., X-^ = XA and X^ = XO. N is 

the number of stations at which the integrand must be 

evaluated. There are (N-l) equal sub-intervals of length H 

in INT = X, T - X n . 
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The accuracy of Simpson's approximation increases with 
decreased H. The pyramided technique progressively computes 
the integral, doubling the number of sub-intervals each time, 
until the value of the computation ceases to change to 
within some desired criterion, EPS. 

Figure 9 shows the change in N and H with each step. 

The coefficients for each step are shown in their relative 
positions according to station. Observe that once a is 
evaluated for a given step, it is used in each succeeding 
step and need not be reevaluated. Further note that the 
first time a Y. is used, its coefficient is 16: thereafter 
it is 14. 

The integrand in (A-III-1) is 

YtjC U ) = U J ~^EXP(CU-B/U)COS(DU+E/U) , (A-III-4) 

A 



and it may be shown that 

Y r ' (U) =U J_ ^EXP(CU-B/U) jcOS(DU+E/U)(J-5/2+CU+B/U) 

( A-III-5 ) 

- (DU-E/U ) SIN(DU+E/U )J . 

The integrand in (A-III-2) is 

T * 

Y (U) = U *EXP(CU-B/U)SIN(DU+E/U) (A-III-6) 



and 
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' t_ I r 

Yj (U) = U ^EXP (CU-B/U ) jSIN( DU+E/U) ( J-5/2+CU+B/U) 

(A-III-7) 

+ (DU-E/U ) COS (DU+E/U )J . 

STEP ANALYSIS 

(Refer to the subroutine at the end of this section.) 
Subroutine ZSUM arguments are as follows : 

XA = Lower limit 
XO = Upper limit 

J = Summation index set in V/SUM. [Refer to Eqns . 

( 11 - 30 ,- 31 ).] 

The following parameters are brought to ZSUM via common 
space HOLD: B,C,D,E,EPS, and LIM. B,C,D, and E are computed 
in MAIN [Refer to Eqns. (A-I-12,-13,-1^ , and -15).]. EPS 
and LIM are set in V/SUM. XA, XO, and J are also set in V/SUM. 
The values of (A-III-1, -2) are returned on SUMR and SUMI . 

Equations (A-III- 1 / ,-5 ,-6 , and -7) are first evaluated 
at the left endpoint. AA = CU-B/U is the argument of the 
exponential function. If AA is less than (-172), it is 
assumed that the exponential term drives the integrand and 
its derivative to zero by branching to statement 6 (Computer 
underflow errors are generated for smaller values of the 
exponent.). BB = DU+E/U, the argument of the SIN and COS 

f 

functions, and CC = DU-E/U, a term used in Y , are computed. 
The evaluations of Y R (XA) and Yj(XA) are obtained on YAR 
and YAI from the arithmetic statement functions YR and YI 
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where YR is (A-III-4) and YI is (A-III-6). EE and FF are 
additional terms used in the derivative. Finally the values 
of Y R (XA) and Y^ (XA) are computed on YARP and YAIP. 

Statement 6 sets all values Y R (XA), Y^CXA), Y R (XA), 

t 

and Yj (XA) to zero for small values of the exponential 
function (those cases where AA <_ -172). 

The same process is executed for the right end point. 

Here Y R (XO) = YOR, Y-^XO) = YOI, Y R ’ (XO) = YORP, and 
Yj (XO) = YOIP . 

The first step in computing the internal points is to 
set the interval of integration in DINT = XO - XA. Next, 
various summation parameters are zeroed including SOR and 
SOI, the sums of the Y R ^ 1 s and Y^'s which will be multiplied 
by the coefficient 1 4 [Refer to Eqn. (A-III-3)] and SER 
and SEI, those sums to be multiplied by 16. SP sets the 
number of subintervals in each step (initially set at 2) 
and K sets the number of terms in SER and SEI to be computed 
in each step (initially set at 1). 

DO LOOP 100 controls the "pyramid" level or step from 
one to LIM. H -is the length of each subinterval in each 
step, H = DINT/SP. 

On each step, the Y^’s computed the previous step for 
coefficient 16 (summed in SER and SEI) are added to those 
with coefficient 1^ (SOR and SOI). SER and SEI are then 
zeroed in preparation for the next computation. XE is the 
dummy variable in this portion of the calculation. The 
initial Y^ is computed at XE = XA + H. 
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DO LOOP 10 computes, the K new Y^ ’ s for SER and SEI. 

AA and BB are computed as before. Y 0 (XE) and Y T (XE) are 

n ± 

computed on YER and YEI and summed on SER and SEI respectively. 
Finally the dummy variable is augmented by 2H and the loop 
continues . 

SUMR and SUMI hold the values of (A-III-1) and (A-III-2) 
for any given step. TEMPR and TEMPI are used as temporary 
storage for the values in SUMR and SUMI computed in the 
previous step. Once a new set of values are computed in 
SUMR and SUMI, they are compared with the previous values. 

If both do not differ from the previous values by more than 
EPS, convergence criteria is met and the loop terminates. 

Should the process reach LIM iterations without convergence, 
a warning is printed and the existing values of the integrals 
are accepted. 
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A-IV. SUBROUTINE COR 



Subroutine COR controls the computation of the correc- 
tion integrals (A-II-5) and (A-II- 6 ) for the case X = 0. 

The integrals are 

AL T _r 

CCOR = / U *COS(DU+E/U)dU (A-IV-1) 

0 

and 

AL j 5 - 

SCOR - / U° * SIN(DU-» E/U)dU (A-IV-2) 

0 

where D and E are specified in ( A-I-14 , -15) . An analytic 
approximation for these integrals exists in theory; however, 
computer numerical error prohibits reliable use of this 
technique for J >_ 5 and R > 2 [Refer to the Discussion of 
Results for an explanation]. For this reason, the approxi- 
mation is used only for J £ 4. This range includes those 
troublesome integrals having integrands with infinite ampli- 
tude near the origin. Standard numerical integration is 
used for J > 4 where the amplitude of the integrands becomes 
small near zero thus enabling one to approach zero sufficiently 
close to ensure that the omitted interval is negligable. 

CASE FOR J < 4 

Under the assumption that, for a properly selected AL, 
as U approaches zero the term E/U dominates DU in the 
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arguments of SIN and COS then equations (III-40) and 
(III-41) are reasonable approximations of (A-IV-1) and 
(A-IV-2) even when a reduced number of terms is included. 

Such a selection of AL is AL - R/100 which ensures 
( E/U )/DU _> 1.0E4 for U < AL. 

Subroutine COR carries out the operation of (III-40) 
and ( 111-4 1) where 

AL t _£ AL t 3 

CCOR = / U * COS(E/U)dU - D / U *SIN(E/U)dU 

0 0 

n 2 AL t i n3 AL t * — 

-- 3 - / U 2 COS(E/U )dU + 4- / U ^SIN(E/U)dU 

^ 0 0 0 

( A-IV-3) 

and 

AL t__s" AL j j 

SCOR = / U d * SIN (E/U)dU + D / U *C0S(E/U)dU 

0 0 

D 2 AL J — - AL tj.J. 

/ IT 2 SIN ( E/U ) dU - - r f IT 2 C0S(E/U)dU. 

^ 0 0 0 

( A-IV-4 ) 

Due to the inaccuracy in the computation of the integrals 
by Subroutine RED, a power of U greater than 1.5 is avoided 
by reducing the number of terms in the truncated series by 
one each time the series is used. That is, for J = 1, four 
terms are used; for J = 2, three terms are used; and so 
forth until only one term is used for J = 4. 
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CASE FOR J > 4 

For this case Subroutine ZSUM is used to numerically 
compute the integral over the interval XA to AL where XA 
is picked as close to the origin as computer time and 
accuracy warrant (XA = 1.0E-4 was seclected for these 
results). It should be noted, however, that in some cases 
the value of the integral over this interval may be less 
than the convergence limit, EPS, used in ZSUM in which case 
the results may suffer from inaccuracy. Provided the error 
is on the order of EPS (EPS should be the upper limit) such 
inaccuracies are not significant and should not affect the 
end result. 

STEP ANALYSIS 

(Refer to the subroutine at the end of this Appendix.) 

The inputs to Subroutine COR include: 

= Upper limit of the correction integral 
= Summation index in equation (III-30) and (III-31) 

2K 3 

X 2 + 4K 2 

KR 2 
2 

= COS(E/AL) 

= SIN (E/ AL) 

AL, 

= / — r - SIN(E/U)dU 

0 U 

AL , 

= f Ar COS(E/U)dU 

0 U** 



AL 

J 

D 

E 

C 

s 

ss 

cc 



126 



D and E are fixed parameters computed in MAIN; C, S, 

SS, and CC are formed in Subroutine WSUM. CCOR and SCOR 
return the values of (A-IV-3) and (A-IV-4) to WSUM. 

The first step is to determine if J = 1, 2 < J <_ 4 , 
or J > For the case J = 1, four integrals are computed 
in (A-IV-3, -*0 of the form 

CINT - / U n * COS (E/U ) dU (A-IV-5) 

0 



and 

AL »«_ f 

SINT « / U ^ SIN(E/U)dU (A-IV-6) 

0 

where M = J, J+l, J+2, and J+3. These integrals are easily 
computed in four step DO LOOP 10 utilizing Subroutine RED. 
RED reduces (A-IV-5, -6) to the general form 

Integral = f (AL,N,E) + ~ ^™(E/U)dU. (A-IV-7) 

0 u^COS 

The values of these integrals are stored in two arrays, 
SI and Cl, of dimension four. Note that should a full four 
term truncated series be desired for J > 1, only one addi- 
tional integral is needed in (A-IV-3, - i O each time. That 
is, in terms of the J of the preceding calculation, values 
of CINT and SINT are required with M = J+l, J+2, J+3, and 
J+^. Since the first three are already on hand, DO LOOP 
13 sets SI(I+1) into SI(I) and CI(I+1) into CI(I). The 
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fourth position is filled with the newly calculated values 
of CINT and SINT for M = J+4. However, assuming error in 
RED for J > , it was elected to truncate the series by 

one additional term for each case, J = 2,3, and k. Hence, 
zero is set into SI (A) and 01(^1). 

One assumes that this reduction in accuracy is still 
acceptable compared to the error introduced by an attempt 
to use numerical integration. 

In the case J > ^, a branch is provided to statement 21 
where subroutine ZSUM is used to calculate the integral 
between AL and XA. Although this procedure simplifies 
the handling of COR results by V/SUM, inaccuracies of the 
order of EPS are expected as previously mentioned. One 
might expect greater accuracy using another technique 
which did not compute the integral piecemeal. 
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SUBROUTINE COR ( AL , D T E , J , C , S t S S ♦ CC , SCOR T CCOR ) 
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CALL ZSUM( XA*X0, J*CC0R.SC0R) 

SC0R=-S COR 

CONTINUE 

RETURN 

END 



A-V. SUBROUTINE RED 



Subroutine RED provides the value of the general integrals 



and 



A t-E 

SINT = / U d *SIN(E/U)dU (A-V-l) 

0 



A J-£ 

CINT = / U z COS (E/U')dU (A-V-2) 

0 



by reducing them to functions of the form 



SINT = F (A,J,E) + / 4 pS(E/U)dU (A-V-3) 

1 o U * 

and 

A 

CINT = F 2 (A,J,E) + / lj™(E/U)dU. (A-V- 4 ) 

The integrals in (A-V-3, -4) are functions of the Fresnel 
integrals as follows: 



CC = / COS(E/U)dU = - C (E/A ) ] (A-V-5) 

0 U* V L 2 

and 

SS = / A ^SIN(E/U)dU = - 5 (E/A)] (A-V-6) 
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where 



E/A 

C(E/A) = / ----- - COS (T) dT (A-V-7) 

0 /2-ttT' 

and 

E/A 

S(E/A) = / SIN (T )dT. (A-V-8) 

0 v^ttT 1 



That is, SINT and CINT for any J may be reduced by 
integration by parts to functions of the basic Fresnel 
integrals . 

Integration by parts provides the following relations 
for (A-V-l). 

If J = 1, 



CINT 



/ 

0 U 



COS(E/U)dU. 



.( A-V-9 ) 



For J = 2 , 



CINT = 2 [ A^COS (E/ A ) - E / -^SIN(E/U )dU] . (A-V-10) 

0 U* 



And for J > 3> 
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CINT 



- <J -- 2 > 1 {/.'* 



(2J-3)! 



fcOSCE/A) f ' > 

L N=1 



( J-2N )~ 



j 

[ r"' 1] r9M -, vnr .. nN-I-^N-1. J-2N-1 2N-1 
- SIN(E/A) E L2(^-2N-1_)J!(-1) 2 A E 



N=1 



( J-2N-1 ) ! 



+ INT 



'] 



(A-V-ll) 



where 



INT 



2 2J -V“ 1 



(- 1 ) 



f— 1 A 1 

l 2 j ; — COS(E/U)dU, 

0 u * 



(-D J < o 

(A-V-12a) 



or 



INT 



2 2 J -V- 






A 



/ SIN(E/U) dU , (-1)’ 



U 



> 0. 

( A-V-12b ) 



Similarly, it may be shown that for (A-V-2) one obtains: 
For J = 1, 



SINT = / \ SIN(E/U)dU, 

0 \J Vz 



(A-V-13) 



if J = 2, 



SINT = 2[A 2 SIN(E/A) + E / C0S(E/U)dU], (A-V-14) 

0 U /2 
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and for J _> 3, 



SINT " ifjilTT {** ■ 



r ^'^roM ™ni( , xN-l 0 i4N-3 A J-2N Tr 2(N-l) 

I SIN (E/A) Z [2 (J-2N) ]!(-!) 2 A E 



N=1 

[f-.l] 



CJ=2NT! 



ron on tm,/ , vN- 1 0 4N-1. J-2N-l tr 2N-l 

X V f 2 ( J-2N-1 ) J ! ( -1 ) 2 A E 

+ costE/A)^ {jrjN-ryr 



+ INT 



J 



(A-V-15) 



where 



INT 



2 2J-3 e J-1 



(- 1 ) 



/ \ SIN(E/U)dU, 
0 U 2 



(-1) J < 0 

( A-V-l6a) 



or 



INT 






1 ) 



/ COS(E/U)dU, 
0 U * 



(-1) J > 0. 

( A-V-l6b ) 



The notation [ARG] indicates the sign of ARG with the 
largest integer less than the absolute value of ARG. The 
use of [ARG ± 0.1] ensures the proper value is selected and 
excludes the possibility of computer round off error. 



STEP ANALYSIS 

(Refer to the subroutine at the end of this Appendix.) 
The arguments of Subroutine RED include: 

A = Upper limit of the correction integral. 
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E 
J 
S 
C 

ss 
cc 

E is set in MAIN and carried in common to WSUM and as 
an argument through COR to RED. A,S,C,SS, and CC are 
formed in Subroutine WSUM and are carried as arguments to 
RED via COR. J is set by COR. The values of (A-V~l,-2) 
are returned to Subroutine COR on SINT and CINT. 

A test is made to determine if J < 3. If it is, the 
routine branches to statement 3.000. 

If J > 3, equations (A-V-ll) and (A-V-15) apply. II 
and 12 are the summation limits where INT(ARG) = [ARG] 
described previously. The primary series DO LOOP computes 
the first summation in (A-V-ll, -15). The secondary series 
DO LOOP computes the second summation. A test is made to 
check the criteria of (A-V-12,-l6) and proper branching is 
initiated. Statement 31 computes ( A-V-12a , -l6a ) for 
(-1)^ <0, and statement 32 computes (A-V-12b ,-l6b) for 
(-1)^ > 0. Finally, a statement 40 combines the terms to 
form (A-V-ll, -15) on SINT and CINT. 



= KR 2 /2 

= Related to the summation index in (111-30,-31) 
= SIN (E/A) 

= COS (E/A) 

A 1 

= / —x. SIN(E/U )dU 

0 U % 

A 1 

= f At. C0S(E/U)dU 

0 u 3/2 
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For the case J < 3>-a test is made to further restrict 
the case to J - 2 or J = 1. For J = 2, (A-V-10,-1^) are 
used. For J = 1, (A-V-9,-13) are computed. 



136 



SUBROUTINE RED(A,E,J,S,C,SS,CC,SINT,CINT) 
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3000 CONTINUE 

WRITE(6»9000)JtSINT » C I NT 
90 00 FORMAT ( 'O’ , 20X , I 3 , 1 OX , 1 PD 1 5 . 4 , 10 X , D1 5 . A ) 
RETURN 
END 



A-VI . SUBROUTINE IFAC 



Function subroutine IFAC returns the value of the 
argument factorialed, e.g., IFAC(5) = 5* 4- 3*2*1. It has 
provisions for IFAC(O) =1. A copy of the subroutine 
follows . 
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APPENDIX B 



STEADY STATE POTENTIAL SOLUTION 



B-I . COLLOCATION TECHNIQUE 

In order to investigate the perturbation velocity 
potential along a cascade blade for the steady state and 
it phase shift in adjacent blades, one may use the following: 



n ( X -i X 2 

4 W X ’ 0) = t=H 1 ; -—EXPC-^I^Uds 

K 0 /ttT (-1 /iwT -i /5P? ,u . 



+ EXP[- -^yDdS 



\firX* N=0 1 ' -1 /x-s 



X N 

P nM / — - — dS 

0N -l J 



(111-47) 



The second and third integrals are of the form 

X oN ? 

I = / — — E XP [ -B/ ( X-S ) ] dS, where B = XR /4 . 

-1 • (B-I-l) 

Letting 



U = B/ (X-S ) , 



(B-I-2 ) 



one may write (B-I-l) as 



I = /B f \ (X-B/U) N e U dU 
B U V * 

X+l 



(B-I-3) 
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or 



CO 

I = /S' X N / %l-B/XU) N e" U dU. 
B U /<E 
X+l 

The binomial expansion provides 



(1»B/XU) N 



? (-1) K N! ( B/XU ) K . 

K=0 IN=1T!K! 



Hence (B-I-4) becomes 



(B-I-4) 



(B-I-5 ) 



t _ /7v N m , (-1) K B K / .°° TT - (K+f) -U JlT T ^ 

I V B X N ’ Q (N-K)*! K ! X K n f U e dU * (B- 1-6 ) 



B 
X+l 



It may be shown that integration by parts provides 



00 




-X 



— X « 



dX = [-, 



M-l ^m-1 ^2 



1 

M-l 





X m_1 



dX, 



M > 0 



(B-I-7) 



where M may be any half odd integer greater than zero, 

e.g., M = 1/2, 3/2, 5/2, etc. Continuation of the integration 

process yields a general relationship 



~ -x 
/ — r— dX 



X 



m 




(M-I-l) ! (-1) 1 1 1 

(M-l) ! (P 2 ) m_I J 



, (-l) m ?VtP ERFC(P) 
(M-l) ! 
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(B-I-8) 



whore 



oo _u 

ERFC(F) = — / dU = — -[ 1 - erf (F) ] 

/tT F 2 U 4 /“t? 



(B-I-9 ) 



For the case N = 0, (B-I-6) is 



I = /If / e 



-U 



-1 U % 
X+l u 



dU. 



(B-I-10) 



The last integral in (111-47) may be written 



X „N N Y N-K, , vK 

; 4^ CS = ^ 



-1 /x^s" 



( B-I-ll ) 



and observe that for the case N = 0, (B-I-ll) becomes 



/ _i_ dS = 2/X+T. (B-I-12 ) 

-1 /T-^S 1 

Combining Eqns. (B-I-6, -10 ,-11, and -12) in (III-47) 
and replacing N with J-l and K with K-l, one arrives at 
the final computer compatible equation 



ip v _ n (X,0) = -[2 /x+T + 2/B 1 SSUMC-3/2, b/(X+l)] 

K_0 

i i 2(-l) K-1 (J-l) !X J_K f T1 (X+l) K-1 /x+r 

~ k Z 1 (K-1) ! ( J-K) ! \>,J-1 2K-1 



- v 1 j _ 1 b k_1 /b‘ SSUM[-(K+|) ,B/(X+1)]}] (B-I-13) 
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where 



SSUM( A,FSQ) = i / X A EXP(-X)dX. (B-I-14) 

F3Q 

SSUM is computed in Subroutine SSUM and will be dealt with 
in Appendix B-I-B. 

B-I-A . THE MAIN PROGRAM 

The purpose of MAIN is the computation of (B-I-13). 

Inputs to the program are the controlling parameters, blade 

spacing, R, and thickness parameter, X, as well as the 

values of the interference coefficients obtained from the 

collocation coefficient solution of Appendix A. 

' \ 

Step Analysis 

(Refer to the routine at the end of Appendix B-I.) 

Double precision is used in Subroutine SSUM and single 
precision is used in MAIN. Double precision arguments to 

v 

and from SSUM are carried on variables beginning with W. 

The inputs DLAM = X, DK = K, and RHO = R, and the 
coefficients of the upper and lower interference functions 
^lm an< ^ m = Ij2,...,10) are read into the routine from 

data cards and printed out immediately for verification. 
(Note that K, the reduced frequency, is used for accounting 
purposes on the output only.) Array DEL(20) carries the 
coefficients with the lower blade coefficients in the first 



ten positions. 



The inputs positioned, various fixed parameters are 
formed and the station, X, is initialized. 

B = AR 2 /4 

BH = /F 
C = 1/ /rrT 

D = -3/2 

X = -0.9 

DO LOOP 100 is the station control loop which generates 
a potential value at each of 20 positions on the blade 

(X = -0.9, -0.8, ... , 1.0). 

In preparation for the use of Subroutine SSUM, the lower 
limit is formed. 



WPSQ = B/(X+1). 
SSUM is called to form 



SSUM[D,B/(X+1) ] 



1 

2 



/ 

B 

X+l 




dU. 



With AA = /X+!T, the first two terms in (B-I-13) are formed: 



BB = 2 y/X+T + 2/B 1 SSUM[-3/2,B/(X+l)]. 



DO LOOP 200 accomplishes the J summation in (B-I-13) 
while LOOP 300 carries out the K summation. (Recall that 



L = 10.) Each term in the summation is built up in a 
straight forward manner. Subroutine SSUM is required to 
compute 



SSUM(A.WFSQ) = | WFS q^ U e dU 



where 



A = -(K+|). 

Function subroutine IFAC is also used. This subroutine 
produces the factorials of its argument and is explained 
in Appendix A-VI . 

The summation loop completed, the summation term is 
combined with BB and C to form POT = of equation 

(B-I-13). Finally POT is printed together with its station. 

\ 

B-I-B . SUBROUTINE SSUM 

The purpose of SSUM is to produce the value of the 
general integral 

00 

\ f X A EXP(-X)dX, A = ±1/2, ±3/2, ±5/2, (B-I-15) 

d 2 
F^ 

Only the derivation of the procedure for A < 0 will be 
presented in detail here. The case for A > 0 is derived in 
a similar manner. This latter case is not of interest in 
the potential solution. 



The on line functiop subroutine, INT(A), generates the 
sign of A with the largest integer less than A, e.g., 

INT(-3/2) = -1. The subroutine, ERFC(F), produces the 
compliment of the error function using F as the lower limit. 
It should be understood that (N/2)! = |- • (^ - 1)* (^-2)-.. i 

and that (-1/2)! = 1. 

Now for demonstration purposes, (B-I-8) is expanded 
for the cases A = -3/2 and A = -5/2. 

A - - 3/2 

1 “ -X i , ? 

~ / ^-rr dX = | -2/nERFC(F)] (B-I-16) 

* p 2 X ^ e ? d * 

A = - 5/2 



1 

2 



» X 
/ ^7- dX 



X 



V 2 





2 2 7 ^ 2 2 v/5rERFC(F) 

37IFJ + 3U ] 

(B-I-17) 



A routine may be generated in the following manner. 
Observe that if one defines I = INT(A), there are I terms 
in the series (Those terms not including the error function 
term.). The numerator is increased by two in each succeeding 
term. The denominator of the first term includes (21-1); 
each succeeding denominator in the series is multiplied in 
turn by (2I-l)-2, (21-1)-^, .... The power of F in the 
denominator begins with (21-1) and is decreased by two in 
each succeeding term. The sign changes with each term. Also 



1^7 



note that the numerical^ coefficient of the error function 
term is the same as that for the last series term with the 
sign reversed. 

Step Analysis 

(Refer to the subroutine at the end of Appendix B-I . ) 

p 

The inputs to the subroutine are A and PSQ = F . The 
value of (B-I-15) is returned on SUM. In the first steps 
various fixed parameters are formed, the summation parameter 
is zeroed, and INT(A) is determined. A test is made to 
determine if A > 0 or A <_ 0. Here only the branch A <_ 0 
is of interest although the structure of the unused portion 
of the subroutine for A > 0 is similar. 

For A <_ 0, statement 10 tests for the special case 
A = i (or I = 0). This case cannot be generated by the 
algorithm used for other cases. If A = - ^ , statements 12 
and 14 set 

SUM = ~ ERFC(F) 

and the routine returns to MAIN. 

If A < -jj, statement 13 begins the algorithm. After 
several preliminary steps, DO LOOP 11 forms the series 
pointed out previously. The method takes advantage of the 
observations made in the introductory material. The series 
formed on SUM, it is divided by the exponential and the 
error function term is added. 



The result is returned to MAIN. 
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M= I + 10 

UIP.ITEl6.15)KtDEL(I)»K,DEL(M) 
16 CONTINUE 
WR I TE ( 6,31) 
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SUBROUTINE S SUM ( A , FS 0 , S UM ) 
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B-II. SONIC WIND T-UNNEL VJALL INTERFERENCE TECHNIQUE 



Development of the Sandeman transonic wind tunnel wall 
interference solution for the steady case resulted in the 
following perturbation potential function from Section III-E: 



i|>(X,0) = - 
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imtA 
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X 



v(S) 
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dS 



+ 2 l [ -p- i - j j- EXP[-A 
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,, 2,2 

N d 



dS 



(III-60) 



To compare this solution with the cascade problem, (d) is 
related to blade spacing R by d = R/2 [Refer to Fig. 7]* 
Further, as shown in Eqn. (IT.I-18), 



v(S) = (X, Y = 0) = -1, for K = 0. 



It follows that the potential may be written: 
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(B-II-1) 



Integrating the first term and making the proper substitution 
in the second term one may write 
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The basic integral within the summation may be expressed in 
terms of the error function, i.e., 



/ hj du = 2_|_- 2 ^ ERPC(P) 

u 



where 



•u>2 



_ AN 2 R 2 

" TTx+I) 



(B-II-3) 



(B-II-4) 



and 

ERPC(P) = 1 - erf (P) . (B-II-5) 



Hence, the potential function: 
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(B-II-6) 



Computation of this function is easily carried out by 
computer. The series is truncated when the inclusion of 
additional terms only changes the sum by a value less than 
the convergence criterion, EPS. 

Step Analysis 

(Refer to the routine at the end of this appendix.) 

The inputs to the routine are DLAM = X and RHO = R. 

DK = K is inputed for accouting purposes in the output only. 
The output includes the above parameters and the value of 



15^ 



the perturbation velocity potential at each of twenty 
stations along the blade. 

The control parameters and various fixed parameters are 
calculated prior to entering the control loops. DO LOOP 100 
provides control to calculate the potential at each of the 
stations beginning at X = -0.9 and increasing by increments 
of AX = 0.2. DO LOOP 50 computes the series summing terms 
until convergence is reached. Again, the function subrou- 
tine ERPC is a standard on line library routine performing 
the operation in (B-II-5). 

Once the series is computed, statement 51 forms the 
potential on POT which is subsequently printed along with 
the applicable station location. 
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